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DETIipjOTATIOlT OP OPTIML CUTTING- PARAMETERS AND 
GROUPING OP OPERATIONS OH MUIT I -SPINDLE AUTOMTICS 

A solution methodology is proposed to determine 
optimal cutting parametei-s on multi-spindle automatic 
lathes with minimization of C 3 rcle time as evaluation cri- 
terion. The difficulty in attacking this prohloin arises 
from the following reasons ; 

1. Information regarding grouping of operations is 
necessar 5 '' for optimization of cutting parameters. 

2. This information is not available initially. 

3 . Besides precedence (technological) constraints 

the information about cutting parameters have also 
to be used for grouping, because operations which 
form a group are forced to assume the same ma,gnitud 
for some of the parameters depending on the cons- 
traints of the machine-tool. 

This interdependence between grouping of operations 
and cutting parameters optimization is sought to be over - 
come using a 3 phase procedure. 
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Phase 1 : Optimization of cutting parameters 


considering the operations to be machined 
sequentiali 3 r. The output of this phase 
provides ah input to the nest phase. 

2. Phase 2; Grouping of operations based on the 

cutting parameters generated in Phase 1 
and the precedence constraints. 

3 . Phase 3; Optimization of cutting parameters taking 

into account the constraints generated by 
the grouping of operations in phase 2 with 
minimization of cycle time as the objective 
function. 

Optimization in pha.ses 1 and 3 is camied out using 
Surrogate Geometric Programming, a relatively new tochniQue, 
which takes advanta,ge of the posynomial formulation of the 
problems. In phase 2 a quantitative approach based on 
Cluster Analysis is used for the grouping of operations. 

The method ologjr has been implemented on IBM 7044- 
Por a Sample problem given by Ba.rtalucci et al using the 
Same input data for the test problem, phase 1 yielded 



20 varia^blss and 41 constrs-ints while phase 3 involved 
9 variables and 29 constraints for a 3 spindle automatic 
lathe. It was observed that the proposed approach yields 
speedier results as compared with the results reported by 
Bartalucci et al in terms of lower production cycle 
time and higher computation efficiency. 



CHAPTER I 


INTRQDUCTIOIT 


1.1 TEB PROBLEM - AN OVERVIEW; 

The constantly increasing costs of lahour, machine 
tools and materials and growing competitive activity, hoth 
domestic and foreign, are placing new challenges before the 
metal cutting industry. Metal forming by the chip producing 
methods must be done, not only faster, but at a lower unit 
cost, if profit is to be realized. Inefficiencies in 
machining practices had been and is being tolerated 
disregarded by many for years, because no real squeeze has 
been felt. The era of such complacency has to end for most, 
and the time will rapidly runout for the few remaining. 

This thesis is concerned with the problem of deter 
mining optimal machining conditions on multi-spindle auto- 
matic lathes. Niven a job to be machined on a multi-spindle 
automatic, it is necessary to determine the optimal machining 
parameters of speed and feed for each of the required 
operations and to assign the operations to the availa 
spindles, the latter being referred to as grouping of 
operations. For better understanding of the problem, the 
important characteristics of multi- spindle automatic 
lathe are described next. 
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A horizontal multispindle automatic lathe has a 
number of stations (or spindles) available for cutting 
operationsj on each of them one axial and one radxal 
carriage are available. Cutting tools are mounted on the 
carriages according to the machining requirements and the 
grouping of the operations on the spindles. 

The workpiece is progressively machined at each 
.t.ticn and i. Indexed to the next etation only after machin- 
ing is finished at each station and the carriages carryrng 

the cutting tools Irate heen withdrawn from the worhpieoe. 

-fm-' 1 oadins and unloading which is 
One station is reserved foi loaaing 

performed either manually or mechanically. 

AS far as the spindle motor is concerned, the following 

schemes are in vogue. 

Single power unit driving, all the spindles at the 
same speed 

Single power unit driving all the spindles at 
different speeds 

Separate power unit for each spindle. 

The most commonly used are schemes 1 and 2, h' 
scheme 3 is used in special machines and in transfer lines. 

coming bach to the problem at hand we notice that 
fo, the optimisation of cutting parameters of the various 
operations one needs infort^tion regarding the grouping of 


1 . 


2 . 


3. 



operations. However, the grouping of operations is to be 

« 

based not only on the precedence requirements of the 
various operations but also the values of the cutting 
parameters selected for various opera,tions . This is due 
to the fact that while constituting a group of operations 
for Si station (spindle) the following important aspects 
must be borne in mind. As far as possible, 

1. the group comprises of operations having the same 
optimal speed (rpm) 

2. for operations to be performed by the axial slide, 
the optimal feed rates should be the same for all 
the operations 

3. for operations to be performed by the radial slide, 
the optiroal feed rates should be the sa,me for all 
operat ions . 

In order to overcome this difficu;^.ty which arises 
due to interdependence of the cutting parameters and the 
grouping of operations, a three plmse procedure has been 
evolved. 

In phase 1 the problem is fomulated assuming that 
all the operations of the given job are machined sequentially 
and the evaluation criteriotn is optimised. The output of 
this optimisation, i.e., the optimal speeds and the feeds 
for individual operations, are then exploited for the 
formation of groups in phase 2, 



In phase 2 the output of phase 1 is used to form 
fea,ture vectors, which describe the characteristics of each 
opera-tion. These feature vectors are then grouped using a 
quantitative technique called Cluster Analj^sis. The output 
of this phase then is a grouping of operations which 
provides the input for phase 5. 

In phase 3 the final optimisation is made after 
reformulating the problem so as to account for the grouping 
of operations. The output of this phase will provide the 
optimal cutting parameters. 

Surroga,te Geometric Programming techniQ.ue is used 
for optimisation in both p'hase 1 and p’hase 3 to take 
advanta.ge of the posynomial structure of the formulation. 

1.2 OVERVIEW OP THESIS; 

The main body of the thesis is concentrated in 
Chapter III and literature survey is given in Chapter II. 

In Chapter III, p^roblem formulation and solution 
niethodology are presented. A relatively recent optimisation 
technique. Surrogate Geometric Programming is presented as 
a solution methodology for optimisation in |Jhase 1 and 
phe-se 3* The core of this algorithm is a geometric program 
suggested by Blau [5]> while the scheme for updating the 
surrogate constants was developed in the course of this 
work. The solution methodology for grouping of operations 



5 


using a q.uanti tatiTe technique called cluster analysis 
is also presented. The txTO solution methodologies were 
computerised and implemented on IBM 7044. 


The proposed methodology was tested on a test 
case used by Bartalucci et al.[3]. The same set of input 
paiameters were used. The test case requires ten operations. 
Therefore, phase 1 foimulation involves twenty variables 
and fortyone constraints. The results of the application 
of the three phase procedure are discussed in Chapter IV. 



CHAPTER II 


LITERATHRS SURVEY 

In this work a brief review of the work reported in 
the literature regarding cutting V3,riables optimization 
8 .nd grouping of operations (on multi-spindle o.utomatics ) 
is presented. Since tool replacement policies are 
closely linked with cutting varia,ble optimization ^ some 
of important contributions in this area are also surveyed. 

2.1 OPTIMIZATIOl'T 

2.1.1 Optimiza.tion of Cutting Variables : 

The problem of selection of optimal machining 
conditions in the basic turning opeiations has been analysed 
until recently with vo^rying degrees of generality by many 
investiga.t ors . Gilbert [28] 5 Brewer and RuedOu [8] studied 
the economics of ma,chining and considered the influence of 
cutting speed, whereas, other authors, Brower [7]i Brown [9] 
Brewer and Rueda [8], Cook [10], analjT-sed the effect of 
both cutting speed and feed r?-ta. The afore-mentioned 
analyses were based on the minimization of unit cost. Unit 
cost and production rate were also optimized, by Eield 
et al [24], by determination of the optimal value of 
cutting speed for turning, milling, drilling, reaming and 
to.pping. 
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Brewer [7], Radford o-nd Bichardson [59j considered 
the prohlem of optimizing unit cost in turning taking into 
account the cutting power available as a constraint and 
Petropoulos [54] studied the constrained problem of the 
selection of optimal machining conditions in face milling 
by using unit cost and production rate as optimal 
criteria. 

G-oranskii [50] » ia a. rigorous mathematical analysis 
expressed the limitations imposed by ma.chino tool, workpiece 
and cutting tool as power transforms, solving the relation- 
ships by linear programming. 

The first use of com'puters in optimization seems to 
hci-vo been reported by Doris et al [15] in which Si computer 
program was used to evaluate all combinations of speed, feed 
and depth of cut, seeking the minimum overall cost. The 
limitations imposed by tool wear, power, and workpiece 
rigidity were o.lso taken into account. 

Brown [9] developed an analytical method for studying 
the effect of various parameters in turning in two passes. 

Brewer [7] suggested the use of lagrang* multipliers 
for the optimization of the constrained problem of unit 
cost, with mainly cutting power as constraint, but did not 
proceed with the analytical optimization procedure. 



Bha.ttaclaa.ry3. et 8.1 [4] optimized the unit cost 
for turning subject to constraints of surface roughnoss 
and cutting power by tho use of Lagreungo's method of 
coupled extremes. This method yielded rO. system of highly 
nonlinear equations the solution of which wa.s effect od 
through a computer,. 

Optimization in practical workshop conditions was 
considered by Ham and Gonzalez [33] by itcratiroly updating 
the tool life equation from actual workshop production 
rosults till convergence is obtained and then tho optimal 
conditions are determined. 

Garina et al [27] proposed experimental determination 
of optimal cutting conditions with the help of Simplex Motho 
(Holder and Mead), 

Cost optima.1 speed for .¥.0. lathes was determined 
by Yulfson and Deryabin [7'4] by first considering the feed 
to be constant e.nd then cost optimal feed was found by an 
exhaustive calculation over a3-l poss.iblo food vc.lues. 
Gurevich [38] later modified this 3.1gorithm to take into 
account the Taylori zed range. 

ileznikov [68] optimized time for a single operation 
with speed and feed as variables using line-ar programming. 
The constraints considered were speed, feed, horse power, 
surface finish and quality of sub-surface layer. 
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Jxops [43] studied the effect of constant cutting 
speed and consto.nt rpm on tlio cost and tino parameters 
when mo.chining stepped parts with, largo and small v.?.ri.ation 
in di'c,moters, while Crookall and Venkataramani [11] 
optimized speed, feed and depth of cut with respect to *'■ 
time for multi... pass turning with power, deflection, 
surface finish raid tool wc.a.r constants. 

To bridge tho gap bo two on thooroticall y optimal raid 
practic.ally optiraal -valuos two f'.pproa.chGS .are av.ailablo, 
namel 3 r, detormination of correction factors exporinonfelljr 
and incorporation of the probabilistic nature of tool life 
and cutting force into tho optimization procedure. 

Thus Kolov [42] suggests o, graphical, nothod to 
dotejrmino the spoed and feed 'bj moans of a cyclogram 
constructed for a particular tool and workpiece matoria.1 
combination, and Romanov and Sotnikova [69] present correc- 
tion factors obtained after extensive experimentation 
considering factors such o.s sot-up, part-clamiping, etc. 

Iwata et al [36] proposed an analsrtical method 
applying a ch.ance constrained .approach to determine cost 
optiiml cutting conditions - speed and food - considering 
tho probabilistic nature of objective f-unction (tool life 
exponents) and constraints (cutting force exponents). A 
solved Gxajnplo indica.ted a, ro-ther insignificant change 
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- 2.594 yen/pieco ( 1 paise/piGCC ) - for a cli'.nge in 

proToability level of tool life from 0.5 to 0.99. 

In. a. more recent work Iwata et al [57] optimize 

the Speed, feed and depth of cut with respect to cost. 

Depth of cut is ■ ;f irst optimized using djrns.mic programming 

after which a chance c ons t ro.ine d approach and SUM! arc 

used to optimize the speed and feed as in the previous [56] 

work. Here again one notices an insignificant cho.nge , tho/t 

too only in the feed, of 0,0053 mn/rev for a variation 

in prohahility level of the constraints from 0.5 to 0.99. 

This is insignif icc-nt in the sense that unless the machine 
of 

as meant to he/high precision, such precision is not normally 
available. Effect of tool life (objective function) 
probability level variations were not presented. 

Geometric Programming (GP) has also been used b3'- 
m,any investigators for the determination of optimal cutting 
parameters. Since the present woirk also uses G P, an 
extensive review for the same is presented in the following 
section. 

2.1.2 Geometric Programming ' 

2.1. 2.1 The GP Techniques . 

The first extensive work on Geometric Programming 
was brought out by Duf fin, Peterson o.nd Zener [17]. In 
the la,tcst work on the subject Phillips and Boightier [58] 
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covGrod all the Torsions of Q-Gomotric Programming avahlahle 
till 1975 a.nd discusoedtho merits and demerits of each 
along with a large number of engineering and other a.ppli- 
c at ions . 

Ragsdell and Phillips [60] used the cond ens.at ion 
version of G P to design o, class of welded structures > 
which problem consisted of 9 nonlinear constraints, 7 vario-- 
bles and 16 DOD^ The results obta,ined are compared with 
those obtained using 8 diff eront , conventional optimization 
techniques. 

Mancini and Piziali [48] developod a strategy based 
on G.P, using normality and orthogonality conditions of the 
duo-1 problem to define candidamco vertices .and then invest i“ 
gating the solution at each vortex for the purpose of 
optina,lly designing helical springs. 

Dinkel a,nd Kochenberger [14] present e. post- op time.! 
sensitivity analysis procedure using P P o.s the basis, 
in which the effect of chcanges in the coefficients are 
related to the design va^riables. 

Hwner and Mehta. [55] wse the Blau algorithm for 
G.P. to optimize computer pei’formance with throughput 
and utilization as criteria. The problem foimula ted had 
2 variables , one constraint and 7 DOD. 


+ Footnote on page 12 
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Nicoletti and Marian! [50] also used the Blau 
algorithm to solve two separate problems in management“ 
structural control in a graded manpower system and advertis 
ing scheduling - the latter consisting of 5 constraints, 

9 variables and 15 degrees of difficulty (DOD)"*'. 

Rao et al [64] criticized the Blau algorithm - 
because it involves a matrix inversion of a high order at 
every stage of iteration - cud the algorithm developed by 
Templeman et al because onehas to resort towards the 
end to crawling to obtain satisfactor 3 )" convergence - and 
presented a new hybrid algorithm based on Liapunov concepts 

2.1. 2, 2 G-.P. and Machining Varis^bles Optimization '; 

Phillips Sind Beightler [57] first applied G.P. to 
optimize unit cost with speed and feed as variables and 
one, horse power, constraint. The m:aximum degree of diffi- 
culty (DOD) involved was 1. Walvekar and Lambert [75] 
Solved a similar problem but with two constraints of horse'' 
power and surface finish using G.P. Here also the maximum 
DOD was 1. 


+ with respect to the primal of G.P. formulation the 
difference between the number of variables (H) plus 
unity, and the number of terms (T) , including the 
constraint terms, is called degrees of difficulty. 
It indicates the difficultsr involved in solving the 
dual of the G.P. program where there are a, system 
of (N+l) linear equations in (T) variables. 
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Ermer [20] solves the mcliining cost optimization 
problem with G.P.^with constraints of feed, surface finish 
and horsG power. For the problems solved, the mascimum 
DOD involved was 2. 

Petropoulos [54] compares G-.P. with conventional 
methods in optimizing the unit cost for the turning problem. 
The numerical examples involve a DOD of 1. 

2.1.5 Multi-Spindle Automatics : 

The first optimiza,tion effort in this area seems to 
be that of Gilman et al [29] where the optimization tech- 
nique used is not available but essentially consisted of 
first selecting the feeds from reference data. The criteria 
considered were machining time, operation cost and tooling 
expenditure. 

Bantalucci et al [5] optimized the cutting conditions 
on multi-siDindle automatics using a method, which borders 
on exhaustive enumoration, whereby the constrained range 
of variation of each variable is first fixed, softer which, 
starting from the maximum value,' it is progi’ossively 
decrea-sed until the rate of cha^ngo of cost becomes zero.. 

Kapustin et al [39] considered the caso of vertical 
multi-spindle automatics. After grouping the operations 
onto the spindles the optimal speed and feed for each group 
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were determined from the available discrete speeds and feeds 
by complete enumeration. Knowing the group "7hich requires 
the ma. 2 imum cutting time, the speed-feed combinations 
for the other groups were a.d justed to equo.lise this maximum 
time, as far as possible, using a.n assignment algorithm with 
minimization of torque a.s criterion. It may be noted here 
tha.t it is not necessary that discrete feeds be assumed, 
which may be inadvisable in mass production environment 
where cams controlling feeds can be made to order. 

Kothari [44] attempted probabilistic optimization 
of cutting variables on multi-spindle automatics but incorrect 
formulation of constraints invalidated both the approach and 
the optimization results. 

2.1,4 Tool Replacement ; 

Optimal tool replacement policies are closely 
inter-linked with optimal cutting parameters since it is 
the latter which control the tool-lives, which in turn 
influence the optimal tool replacement policy. In this 
section a brief survey of tool-life models and tool replace- 
ment strategies is presented. 

2. 1.4.1 Tool-Life Models s 

The use of probabilistic tool-life models is a 
relatively recent phenomenon the latest of which, in 
available literature, are that of Eamalingam and ¥atson [63] 



15 


who developed a 'single injuip/ tool life model' purely on 
the theoretical basis of the relevant physical wear mechanism 
involved. Hamalingam [62], in continuation of the previous 
work, developed models for 'multiple— in juiy tool li-® » 
again on a purely theoretical basis, but in a diocussion 
Cohen and Black [62] validate EamaJ_ingara' s model by 
comparison >Jith eleven experimental teot repoit.-' 
life. 

A rather interesting model of the mean life of 
equipment subject to compound failure was developed theore- 
tically by Malik [46] who draws the attention of researcners 
toward developing a relevant integral transform of this 
function. 

2.1.4. 2 Tool Replacement Strategies'; 

Kothari [44] mentions four reported attempts at 
o-btaining an optimal tool replacement policy. The first, 
suggestea by Taha, Is a cost model for a strict inter^ral - 
replacement policy where the cost of rejects is taken into 
account. In the second, Ebert and Hershayer obtained an 
optimal with respect to tool - operating costs through 
simulation, the parameters were quality control decision 
rules and strict-interval replacement. Strict Interval 
replacement was reportedly found better. In the third, 
Okushlma, and Eujll formulated cost models for scheduled and 
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unscheduled tool replacement policies and concluded, after 
simulation, that scheduled changes are more economical. 

The fourth by Duncan is connected with mathematical and 
numerical models for scheduling cutting-tool changes. 

Ssterzon and Baclzievskii [21] introduced the concept 
of optimal group-tool changes for which the stoppages in 
transfer lines we re a minimum. This work is based on 
time studies of the tools in a transfer line . 

Sheparov [71] suggested a simple graphical method 
for planning optimal group tool changes which consisted 
essentially of grouping the tools with similar tool-life 
magnitudes into a predet errolned number of groups. 

Pas’ Iio [ 51 , 52 , 53 ] in three separate works, developed 
mathematical models to determine both individual and group 
optimal tool replacement policies. A V/eibul tool-life 
distribution was assumed. 

2. 1.4. 3 Related Topics; 

A tool edge can be used for a number of different 
operations depending on its life. So, when machining a 
number of pieces consequentively on axitomatic machines the 
necessity of optimal tool selection arises. This topic has 
been dealt with by Etin et al [22]. 
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De Vor et al [12] conducted experiments to determine 
the nature of tool-life variation as a function of pre- 
specified wear level and introduced the concept of optimum 
wear level specification as a result. 

2.2 GROUPING OP OPERATIONS: 

Once optimization is completed it becomes necessary 
to group the operations such that they are distributed 
among the available spindles according to some accepted 
criteria. 

2.2.1 Multi-Spindle Automatics: 

Onl37’ three ca,ses of attempts at such a grouping for 
multi-spindle automa.tics were reported in available litera- 
ture. 

Thus Jones and Morgan [38] used the AIDA (Analysis 
of Interconnected Decision Areas) technique to eliminate 
the technologically impossible combinations of operations 
and then examined all the remaining alternatives to determine 
the cost optimum combination, while Bs„rtalucci et al [3] used 
a computer program to enumera.te all the combinations, 
eliminate the technologically infeasible ones and scan the 
rest for the cost optimal one. 

Kapustin et al [39] used a 6 digit binary code to 
progressively group opeiations, the only criteria for 
grouping being that the combination should be technologically 
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2.2.2 C l uster An alysis ; 

This technique of grouping 'objects' which are 
'similar' with respect to some criteria appears appealing 
as it may he possible to suitably incorporate the grouping 
optimization criterion into the feature vector of each 
' object ' . 

Duran and Odell [18] provided an excellent survey 
of different clustering approaches and distance f ujictions , 

while Ha,rtigan [34] and Anderberg [l] made available to 
prospective users different versions of each technique 
along with their computer progrs^ms. 

Duda and Hart [16] discussed the use of various 
clustering algorithms with respect to pattern cla,ssification. 

An important report of computational experience was 
given by Fromm and Horthouse [26] who oAso presented 
empirical relations for lumiping and splitting of clusters. 

A review of cluster aneAysis is given by Kennedy in [41] 

The next chapter presents the problem formulation 
and solution methodology for the problem of optimization 
of cutting variables and grouping of operation on multi- 
spindle automatics. 



CHAPTER III 


PROBLEM PORMUIATIOH AND SOLUTION METHODOLOGY 


3.1 INTRODUCTION 

The type of automatic lathe considered in this thesis 
is as follows 

i) A single motor drives all the spindles. Continuously 
varying speeds are available on all the spindles 
within a specified range. 

ii) Different continuously varying feeds are available 
for axial and radial slides on each spindle. Further , 
the range of feed variation for both types of slides 
are equal. 

The minimization of cycle time is considered as the 
criteria of optimization. 

Given a job to be machined on an automat, it is 
necessary to know the grouping of operations onto the 
available spindles before we determine the optimal cutting 
parameters for the various operations to be performed on 
the automat. However, grouping of operations cannot be 
done without knowledge of optimal cutting parameters for 
the various operations. This being necessitated by the 
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fact that after grouping, in each group, all operations 
have a single speed, all radia,l operations have a single 
feed and all azial operations have a single feed. The 
difficulty caused hy this interaction between grouping of 
operotions and cutting parameters is overcome through the 
use of a 3 phase procedure. An overview of the 3 phase 
procedure is presented nezt. The complete details are 
discussed in Sections 3*2, 3.3 and 3.4. 

Phase 1; In this phase the problem is formulated assuming 
that the operations required for the job are 
ezecuted sequentially and the optimization 
criterion adopted is minimization of total machining 
time. The , following constraints are considered; 

i) iaazimum and minimum feed limits for each 
operation according to the cutting force, 
surface finish s.nd machine tool^ 

ii) mazimum and minimum speed limits for each 
operation a.ccording to the taylorized range 
and surface finish, 

iii) mazimum and minimum spindle speeds of the 
machine tool, and 

iv) tota,l power e.va.ilable on the multi-spindle 
automatic. The tota,l power requirement for all 
operations should be less than or equal to the 
power available. 
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As would be explained in Section 3.2, the problem is 
posynomial in character a,nd therefore Surrogate Geometric 
Programming (SGP) has been used for solving the optimization 
problem. 

Phase 2 ; The optimal cutting parameters obtained in 

Phase 1 are used to construct a feature vector 
which describes the characteristics of each of 
the operations and bluster analysis is used to 
find a grouping of operations talking into account 
the precedence constraints. The output of this 
phase yields grouping of operations. 

Phase 3» The output of phase 2 provides the grouping of 
operations required for the final optimization 
which is carried out in this phase. The problem 
is formulated considering the operations to be 
machined in groups. The optimization is carried 
out considering the va.rious groups of operations 
such that the cycle time is minimized. It needs 
to be pointed out that for the multi-spindle 
automa-tic lathe under consideiution , the various 
groups of operations are machined simultaneously. 
Of the various operations, the operation requiring 
the longest processing time (PT) dictates the 
cycle time (CT) for the job. The processing time 
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is defined as sum of the in 3 .cliiiain. 2 : time and 
the average tool change time per 30 b, for the 
tool concerned. 

The following constra.ints are considered? 

1) All the operations in a particula,r group must have ±h 
the same speed, further, the selected speed should 
be within the specified maximum and minimum speed 
limits of the machine tool. 

2) The feed for all radial operations of a group should 
be equal subject to maximum and minimum limits for 
feed. 

3) The feed for all axial operations of a group should 
be equal subject to maximum and minimum limits for 
feed. 

4) The total power requirement should be less than or 
equal to the total power available on the machine 
tool. 

5) For no opero-tion the processing time should be more 
than the cycle time (CT). 

3.2 PHASE 1? OPTIMIZATION BEFORE GROUPIIG; 

The problem for this phase is formulated assuming 

the operations are machined sequentially subject to the 

various constraints. 
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3*2,1 Objective Function s 


The evaluation criterion is the minimization of 
total ma,chining time which is given by 


"D 


MT = y (t . + t . ) 

. — r- 'mi Cl T. 

1=1 1 


(3.1) 


whore t ^^ ■ = machining time for i-th operation (min.) 


^^i 

S. 

1 


P 

t 


Cl 




^i 

= 

■ 11 

= Length of i-th operation (mm) 

= cutting speed for the i-th operation 
= feed for the i-th operation (mm/rev.) 
= total number of operations 
= tool change time for i-th tool (min. ) 
= life of the i-th tool (min. ) 


The tool life T^ is obta.ined from the following 
expression [3] 


Si d. 


Pi P2 




T. 




'1 


where = a constant, 

d^ = depth of cut for i-th operation, 

Pl,P 2 >P’^aJ® exponents for d^^, f^^ and T^ respectively. 

If eqn. (3.1) is expressed in terms of speed and 
feed wc have 

MT = + Oji S.e f.l^) (3.2) 

■ ' 1=1 - 
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where C^. = 


'2i = \ 


Cp. = (L d, 


) 


a, b, g and h are exponents given by 

1 


a = -1: 


b = -1 


« rr — 

9 O — 




- 1 j h 




The coefficients end C2^ are positive since all 
the terras' in their expressions s-ro positive hence equation 
(3.2) is a posynoraial. 

3.2.2 Constraints; 


The following notation is used; 

ed for the i-th operf 

ma,xiniuin speed limit for the i-th operation 


= speed for the i-th operation, rpn, 


S, . 

tx max 


s 

ti min 


S. 


m max 


S 

m min 
"'i 

f , 

ti max 


according to taylorized range"*” and surface 
finish, rpm. 

= minimum speed limit for the i-th opera,tion 
according to ta.ylorized range and surface 
finish, rpm, 

= maximuTii speed limit of machine tool, rpm, 

= minimum speed limit of machine tool, rpm, 

= feed for i-th opera.tion, ram/rev. 

= maximum feed limit for i-th operation 
according to taylorized ra.nge and surface 
finish, mm/rev. 


+ Taylorized range is the. range for which the Taylor's 
equation is valid. 
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f,., = maxinun speed Unit for nacliino tool.,, 

Licl Ji. 

nn/rov, 

niii = niniriun speed liriit for n.aclaine tool, 
nin/rcv. 

file various cons traiiits e.ro represented as follows* 
1) rfexinun speed constro.int 



(3.3) 


(3.4) 


(3.5) 


( 3 . S ) 




(3.7) 


wliere = power requirement for i-th operation (HP ) 

_ naximum power available on machine tool (HP) 

ricix ■ ' 
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The force required for eo.ch operation [3] is given toy. 



"r 




Ur 


f . 
1 


where K 2 is a constant and paid are exponents of d^ 
and respectively. 


Therefore , 

Ti: D- S . p. 

Hi = -tIooW-^' 


(3.8) 


where 

'0 


diaEietor of tho i-th opeiaation, mn 
efficiency of transmission of povrer to the 


Eg,. 


spindles 


to 


(3.8) is alwo.ys positive with respect /the 


variables and fj^, therefore constraint (3.7) is a 
posjrnoiaial. 


Constraints (3.3) to (3.6) arc- also possrnomials with 
coefficients of unity and exponents of unity/thus, o-ll tho 
const^raints are posynonials . 


3.3 PHiSB 2; IROUPIFG OP OPBRATIOlfS ; 

3 . 3.1 A n Overview ; 

Since one of the objectives of cluster analysis is 
to obtain an insight into the underlying structure of the 
data or to find out * natural groups ' present in the data 
available, cluster analysis has been used for the operations 


grouping problem. 
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The obvious choice of characteristics to describe 
each operation are all data loertaining to speed and feed 
since these two are the control va.riablos. Precedence 
(technological) constraints could be se^tisfied by including 
separo-te character is tics 5 , for each opera/bionj, which would 5 
in effect, keep two operations, which cs.nnot be machined 
simultanieously , sepa?rate. 

Further, the study has been made for 2 cases where 
the clmractcristics of the operations considered a,rG : 

1 ) ma.ximum and minimum speed and feed liBiits, and 
precedence constraints only, and 

2) optimal (Phase 1) speeds a,nd feeds, and precedence 
constraints onl 3 ?'. 

If the first ca.se (or both) provides 'good’ solutions, 

we ca.n group the operations prior to the optimiz?/Gion of 
cutting variables. However, if only the second case produces 
' good ' solutions then the optimization procedure of ihaso 1 
has to precede grouping. 

In the nozt section, the technique of cluster analysis 
is presented, while in section 3.3.3 the details of the 
procedure adopted for the grouping of operations are 
presented. 
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^ ^ ^ __Ana,lys is Teclniiqu3 : 

3. 3- 2.1 Introductions 

In social, economic, industrial, psycliologico.l , 
tiiologica.1 and many other fields it is necessary to group 
similar items together, either to obtain an insight into 
the underlying structure of the data or to classify the 
items under study. The process of forming groups of siriilo.r 
items is colled Clustei’ A nalysis . Thus, given a date, set, 
the objective of cluster anal 3 "sis is to find out 'the na.tural 
groups' present in the data. 'The solutio)! of a clustering 
problera should give the number of clusters present in the 
data and should list the members of each cluster. Obviousljr 
the members of one cluster should, exhibit similarity with 
respect to a set of characteristics or features. 

3.3. 2. 2 Ma.themr,tica-1 Formulation'; 

Suppose tho.t we have a set 0 = [®j .] = 3, , 2 , . . , , M] 
of M objects and another sot C = [C^. | 1 = 1,2,..., h] of 
N characteristics (features) such tho.t eo.ch characteristic 
(C.) is a measurable feature? of eoich object (O-). 

let X- • be the vaJme of the i-th chs.ract eristic (Cj_) 

1 j 

of the j-th object (O^-), l®t Xj = denote the 

I X 1 vector of such aieasuremeiits,. Viewed geometrically, 
the II objects can be thought of as points in an F-dimensional 
space. 
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Given M objects and the data, set X = [X-], a clustering 

J 

algorithm will determine the na.tural clusters present in the 
de.ta.; Such that objects . be3.onging to the same cluster are 
'highly simila.r' . Some 'criterion function' is used to 
evaluate the partitioning of objects. 

3.3. 2,3 Similarity ne8.sure; 

''■'hen do "we say that the two objects i and j are 
similar? how, since the H objects to be clustered can be 
represented as i'l points in an IT-dimensional space, Euclidean 
distanice is the most obvious (and also vridcly used) similarity 
measure. ¥e keep two objects i and j in the same cluster if 
the distance between X. and X. is smaller than some specified 

J- J 

threshold (D^) and in different clusters if the distance is 
more tlia.n D^. If the number of clusters is specified a 
priori then each object is assigned to the cluster it is 
closest to. 


The Euclidean disto-nce between any two points i and j 
in an I'T-dimensional space can be computed from the relation-- 
ship ’ 

I'T o 1 

(3.9) 


d 


13 


( 1 . 1 - 


) 2 1 i- 

^jk^ -I 


where x.^ and x.i are the values of the k-th charact e3?istic 


of the i-th and j-th object respectively, and d^j is the 
dista,nce between them. 
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The difficult]/ in using Euclidean distance as a 
measure of similarity is tlia.t it is affected "by change of 


sca.le. To 


variables 


overcome this drs-wbaclc; we c 
( clmract eristics ) b]'- dividing 


an standardise the 
the k-th characteristic 


by the ma.rimum of all the k-th charax t oris tics j i.e., normali- 
zation. The distance between two objects i 8.nd j then 


b ecomes 


H 




[n ( 

k=l 


^ik 


max z 


“ -iv 2 ± 

]a 


-ttM 


(k 


3.3. 2.4 Criterion Functions 


( 3 . 10 ) 


In a, clustering problem we want to partition a set 
Z = [Z^j ^ 2 , Zj^] of M objects into exactly Z clusters ^ 

where IC ^ Mj, such that the objects belonging to the same 
cluster are similar (close to each other) and the objects 
belonging' to different clusters are dissimila.r. Now the 
number of ways of partitioning M objects into Z distinct, 
nonempty clusters is given by Stirling's number of second 
kind [18 ] 5 

IC 

W = jk st (?) (-1)^ ('K-J)" (3.11) 

U. J 

(e.g,, the total number of ways of partitioning 50 objects 

49 

into 3 dust el’s is appro xima 'be I])- equal to 2 ^). Out of these 
N possible partitions we want to find a partition which 
optimizes a given objective function. The most commonly 
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used objective function in cluster ans.lysis is the sum-of-- 
squared-orror criterion. 

Let be the number of objects in i~th c3-uster and 
be the mean vector of the objects belonging to the i--th 
cluste3r, so that 

Y. = jh ^ ’ i = l,2,...,K (3.12) 

i lUi 


The mean vectors Y^j ^ 2 ^ Yj^ are called 'centroids' 

and for b. given cluster i, Y^^ is the best representative of 
the objects belonging to i. 

ITow the sum--of-.sq,U8,red error in the i-th cluster is 


given by 5 


1 


X 

J 


^i 


(3.13) 


and the total sum-of-squared-eiu’or in all the IC clusters is 


1C 


ii 


s c 3, = z: ^ 


i=l 


i=l -j^" 


X. - Y. 
3 


(3.14) 


This criterion function has a simple interpretation. 
The suia-of-squarod-error (E) measures the total squared 
erroi’ incurred in representing the M objects I 29 ...? X^,,- 

by X cluster centroids xY-j^s Y 25 ... ? Yj^ [16J. 
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5*5.3 Grouping of Oper a tions Through Cluster Analysis s 

The cha,ract eristics used to describe each operation 
aro the maximum o,nd minimum limits of feed and speed, and 
the optimal speed a.nd feed selected in Pho.se 1. This choice 
is dicto-ted by the fact that speed and feed aro the control 
variables. The precedence or technological constraints 
which have to be represented as a. characteristic are discussed 
next. 

Consider the workpiece shown in Fig. 1, this job 
requires 10 turning opera-tions (la.bolcd from 1 to 10 in the 
figure). Of those 10 operations, some cannot be machined 
before others, viz., machining of opcro.tion 2 has to 
precede that of operation 3, operation 6 has to precede 7? 
and operation 9 tmis to precede 10. This co.n be represented 
diagrammatic ally as 



whore the numbers in the circles represent the operation 
numbers and the links represent precedence constro.ints , the 
arrowheads indicating the direction of precedence. Thus, 
for this job there is only one level of precedence, i.o., no 
opcra,tion is preceded by more than one operation. Further 
notice that the unlinked operations 1,4,5 and 8 can be grouped 
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X'j’ith ea.cii other and any of the other operations, for instance, 

operation 1 C0.n be grouped with operation 2 or operation 3. 

that 

The similarity measure used is Euclidean distance so/ all 
these above mentioned precedence features xrere represented 
by introducing an extra chars.cteristic in the fea^ture vector, 
called precedence value, such that if for operation 2 this 
value is 1 then the corresponding value for opera.tion 3 will 
be 0, xhaile that for operation 1 will be 0.5 bcca.use it 
should be equi-distant from 2 and 3. 

The arrangement of the clraracteristics in the 
feature vector representing each operation is 






) 


Xifhere = maximum feed rate (mm/rev) 

0 2 = mimim'um feed fra-te (mm/rov) 

= maximum speed (rpm) 


"4 

‘"6 


= minimum speed (rpm) 

= optimum speed (rpm) 

= precedence value (0 or 1 or 0.5) 


c-rj = optimuEi feed (mra/rov) 


The elements of the feature vectors of all operations 
are presented in Table 4. As indicated in Section p.o.2.o all 
the characteristics are normalized? but the precodence value 
is varied until the groups formed satisfy the precedence 
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constro^int s. This is dono hy multiplying the 'precedence 
values of all operate ions with a, multiplication factor such 
that if 'bhe multiplication fa.ctor is 2, then, for example, 
the precedence values of operations 2, 1 and 3 will hecomo 
2,0, 1.0 and 0.0 respectively instead of the original values 
of 1.0, 0.5 and 0,0 resp ectivol3?-. The value of the multipli- 
cation factor which .just satisfies all the precedence cons- 
traints is determined by tria3. and error. 

Once the feature vectors are formed, groups (clusters) 
are formed using the clustering algorithm given by Mac Queen [l] 
This algorithm is lorown as the IC-me8.ns method [l] - nearest 
centroid sorting - where the nujaber of groups K is specified 
a priori. In our study, K represents the number o'f spindles 
available for machining. The steps in the algorithm are 
as follows'. 

1) Take the first K fea;bure vectors (opora.t ions ) in the 
data set as groups of one menbor each. 

2) Assign each of the remaining M-S, where M is the 
total number of feature vectors (operations), to 
the group with the nearest centroid. After each 
assignment, recompute the centroid of the gaining 


group 
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3 ) After all the feature vectors (operations) have been 
assigned in step 2, take the existing group centroids 
as fixed seed points 3.nd make one more pass through 
the data set e.ssigning- each feature vector (operation) 
to the nearest seed point. 

4 ) hepee.t step 3 till convergence is achieved. (Three 
iterations were found to be sufficient for the test 
problem. ) 

Each grouping opera.tion consists of three separate 
grouping operations as follows 3 

1 ) G-rouping of the radia.l opera-tions using the sx^eed, 
feed' and precedence elements. 

2 ) G-rouping of the axial operations using the speed, 
food and precedence elements. 

3 ) G-rouping of the radial a,nd axial groups found in 
steps (l) e;.nd (2) using onljr the speed and precedence 
elements. Moreover, the elements for this stop are 
the centroids of the groups found in stops (l) and (2). 

The reason for eliminating the feed elements in 
step (3) is that speed is tho 01113;^ common factor, other than 
precedence, betwo on radial and axial operations as separate 
feeds are available for them on each spindle. 

The EORTEA.I'T lY cod a for the K-means algorithm given 
in [34] was modified to incorporate the Mac Queen algorithm 
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and tliG above procedure. The complete listing of the prograeii 
is given in Appendix C. 

In the next section the phase 3 which consists of 
optimization after grouping is presented. 

3.4 PHASE 3 : OPTIMIZATIOH AE'TSH GROUPIEG i 


The grouping of opera/fcions obtained in Phase 2 is 
used for the determination of optimal cutting parameters 
fox* individual groups of operations. The objective 
function and the constraints considered are listed in the 
following sections. 


3.4.1 Ob j ec t lye Func tion : 


The objective function is the minimis, ation of cycle 
time (CT), which is dictated by the opere,tion with the 
maximum processing time (PT). The PT of an operation is defined 
as the sum of the machining time for that opera-tion and the 
average tool-change time per job for the tool concerned. 


The latte 


r u 


3rm having been incl.uded to keep a chock on the 


number of stoppages of the machine tool caused by tool 
failures. Thus the objective fxinction cpai be expressed as 
minimize CT = y 


where y - Max + 


cij 


. mil ') 
T. ^ 


^■*^mik '^ei k T 


^mik 


ik 


)] 


(3.15) 
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3 — 1525,0.5 *^a,i 

1 - _ 1 O p. 

-- ~ -1- 5 5 ... 5 4.0.T 

j. j~ 


All terms inEqs.(3.15) have the same connotation 
as in GO. ( 3 . 1 ) 0 x 0 opt that 

i indicates tho spindle or group number, 

0 indicates the axial opoiation number, and 
k indicates the radia.D. opera/'cion number 


Further, 

h = total number of groups 
^ai ~ number of axial operations in the i-th group 
= number of radial opers-t ions in the i-th group 

Since tho terms in smal3. parenthesis in Eqs.(3.15) 
8 .rG the Same as for Eqs.(3..l) which have already boon shown 
to havG posynoniial structuro, tho objective function in 
this phase also has a posynomial structuro. 


3.4.2 
Let , 


C ons traints, 1 


O 

Uj 



lEc 


speed for the i-th group (or spindle) 

= ma,ximu!n speed for th^ j~th axial operation 
on i-th spiled le 


S. . . = minimum speed for the j-th axial opomtion 

13 mxn 

on ith spindle 
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>0 . 


IK max 


S., 

ik min 


ri 


’ik max 


f 

ik min 


f. . 

13 iin ix 


f. . • = 

13 min 


maximum speed for the k-th Dradial 

operation on i-th spindle 

minimum speed for the k-th radia 3 . 

operation on i-th spixadle 

feed for the re-dial slide on the i-th 

sp indie 

maximum feed for the k-th radial opera tion 
011 the i-th spindle 

minimum feed for the k-th radial operation 
on the i-th spindle 

feed for the axial slide on the i-th spindl 
maximum feed for the j-tli axial operation 
on the i“th spindle 

minimum food for the 3-th axial operation 
on ' the i-th snindlc 


All speeds and feeds are in rpm and mn/rovolut ion 
'espoctively. 

The various constraints are represented as follows; 


1) Maximum, speed constraint 


S. < S 

1 ^ 


where 


i max 
= lan (S. 


a s') 

i3 max^ ik max^ m max^ 

i = 1 , . . . , N 

3 = 1 qg_ 

k = 1 5 ... s ^ y X 


( 3 . 16 ) 
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2) Minirauin speed constraint 


Si > 


s . 

1 min 


(3.17) 


wlior; 


s. .... = ife.x (g, . S.,. 


1 niii:. 


S 


nin " m min 


1 — 1 O « O 0 


j _ 1 , . . . , 


Ic — 1 ^ ^ Q. 


ri 


3) Maninuni radial feed constraint 


f • c< f . 
ri ri max 


where f 


rx max 


Min max’ "^m max^ 


(3.18) 


i — 

j£ =: Q[ 


ri 


4) ilininmm radial feed consura.int 


■^ri ^ ri min 


where - Max min’ ^ min^ 


(3.19) 


i = 3- , . . . , K 

IC = lj,..y l-p-j 


5) I'lo.ximnni axial feed constraint 


f . =$C f ■ 

ai ^ ax max 


whore 


kin (fj_j max’ ■'m max^ 


(3.20) 


i = 1 , . . . , N 

j = 1 , , . . , 9.g_ j_ 


6) Mj.niro.um axial feed constraint 


f . !> f • 
ax ax mxn 


whore = Max (fiJIain^ min 


■i yi ) 


( 3 . 21 ) 


i — . 1 y m ■ f , ■ 

J = 1 , . . . , q.g_ j_ 
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7) Power constraint 




q . 
■^ai 


(Z1 H,. 


i=l 3=1 


13 


q . 

-n 

k=l 


) C H... 


flic 


(3.22) 


x^horo = powor requircaneiit for tiio j-th aicial 

operation on i“tli spindle 
Hij^ = power requiroraent for tiie k-tli radial 
operation on i-tli spindle 

COh-O expression for powo:o is the sarae as given hy 
Eqs.(3. 8) . 

8) llio cjrc3.e time constraints 


for axial operant ions 
*^'^ni3 "'eij 


for radial operations 


4^1 ) 

"i3 


■y 


i = 1, 
3 = 


(3.23 


TT 


' '^-ai 


^ '"mik ‘ 


mix 


cii 




" ) Cy 


Eqs. 


IL -L y o 

k = 1 


(3.24 


TT 


? • • • 9 


■ — 'X— 

Tlio constraints frora/(3.l6) to (5.24) havo posynoniieJ 
structure. The variables invoj.ved are cither speed or feed, 
or both as in Eq.s. (5,22) to (3,24). 
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3.5 SOLUTION HETHODOIrOGy^ 

The structure of the prohloms considcrec: j.n pha.se 1 
a.iid phase 3 are pos 5 ?noinialj, theroforo G-oonictric Progz-auiEiing 
(G-P) has heon used. 

The Blau 3.1gorithm [5] for GP was coupled with the 
Surrogate Geometric Programming (SGP) tcchniqu.G [58]; 
hecauSG the Blau algorithm if used alone, consumes a large 
amount of tine for computation owing to the high order of 
matrix inversion involved. For exa.nple, the oiotinization 
of cutting ps-ranetors for the prohlon shown in Pig. 1 involves 
inversion of matrices of order 62 and 39 for phases 1 8.nd 3 
respectivel 3 E if the Blau algorithm is used a.s such. However, 
when coupled with SGP the order of the matrix to he inverted 
is reduced to 22 in pliase 1 and to 11 in phase 3. 

In the next two sections, i.e. Sections 3.6 and 
3.7, the Blau algorithm and SGP tjchninue will he presented , 

3.6 TEE BMP ALGORITHM; 

3,6,1 The Gene ral Format : 

The most general primal form of a geometric program 
is as follows ; y 
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Minimize y^(X) 


subject to 


= ^n± 


JZ—I 'ot O'C ' i, 
0=1 n=l 


T.. 


m 

^ O’ , C , 
no ni'c 


O'er 


a. 

i;.i oi 


( 3 , 25 ) 






> 0 ; 


m = 1,2 

n = 1,2, IT 


94 -^* « 0 9 


TT 


(3.26) 

(3.27) 


wheri 


^mt = ±. '^ot = -- ^ ^ 

‘^mt ^ ^ot^ ^ 


^'mtn’ "''otii unrestricted in sign 

= number of toms in the ni-th constraint 
(m = 1,2,.,., M) 


= nmber of terms i.n the objective function 


111 terms of engineering design f emulations, the 
C,.,.u are economic coefficients, the X. a.rG design decision 

lu U ii 

variables, the are technological exponents of the 

decision variables, and the a vector has as elenents 


binary variables (+ l), whose signs represent the sign of 
oa-ch term and inequality in the problem stateneiit. 

3 . 6.2 

'Ihe interested reader is referred to [ 5 ] for the^ 
detailed theory, but an explanatory extract is reproduced 
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here ; ’The raain idea is to use a Mcwton-Eaphson proccduro 
[Refs l/ilde and JBeightlor, • 'Foundations of OptiniZrotion' ] 
to orivo to zero the components of tho gra.diont of a Lagrangian 
function formed from the loga,rithms of tho original ohjectiTO 
function and the constraints, A non-linear transforma-tion, 
which amounts to substituting a weighting vsn'^iable for each 
term, makes the Lagrangiann gro-dient lineoa’ in the weights 
as woll 8.S in the La.grangian multipliers, although bilinear 
in the two sets of variables ta.kon together. The poljmonia-l 
form of all the functions mslies available, in closed form, 
the derivatives needed for the Newton-Raphson iteration. To 
begin, one selects a trial priBial solution, not necessarily 
satisfying the constraints exactly. This inniediatoly gives 
tho initial values of all tho weight variables, which, when 
substituted into the Lagrangian gradient, gives a sot of 
functions linear in tho unlmown Lagrangian multipliers. The 
initia.l va-lues of multipliers are then chosen to minimize 
the sm of sq_uares of these functions. Substitution of those 
values of weights and Lagrange multipliers into simple 
formulas gives the numbers needed for a howton-Baphson 
adjustment, obtained by solving a non-singular set of 
linear equations. This gives new values of the Lagrange 
multipliers and of the original primal variables. When the 
latter arc transformed into new weights, the ilowton-Eaphson 
procedure is repeated'. 

The algorithm is reproduced in Appendix 1 
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3.6.3 Users' Ro p o rt s_ ; 

Two reports Imve iDeen recorded in litornturo on the 
use of GOMTRY computer code. In the first, Hwner and 
ilchta [ 35 ] develop a quciiig model for a Computer l/O 
and C.P.U. overlap system and formulate two porfornance 
measures based on tliruput and system utilization. Those 
measures were optimized and validated by appl 3 ?'ing thom t'' 
an existing system. The problem fornulatod had 2 vteriab 
one constraint arid 10 toms. 


In the second report, ilicolotti and Ilariani [50] t 
the computer code to obtain numerical solutions to some 
general problems in management, modeled by discrete nonlir 
systems with nonlinear perforrance indicc^s. Two examples 


disc us s e d relate to 


’structural ccntrol in a graded nanpow. 


system' and to 


' a d ve rt is ing s cb.c d uling ' , 


the latter comprii 


of 9 variables, 23 terms and 5 constraints. 

further, the code was tested satisfactorily on txfo 
of the test problems provided by Dembo [13 ] 5 namely, Heat 
exchanger design problem (8 variable.;, 6 constraints and 
IS terms) and the 3 -stage membrane separation process 
(13 variables, 13 constraints and 53 terms). 
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3.7 SURROGATE GBOICBTRIC PROGRAM s 
5*7.1 Introduction ^ 


A surrogate constraint is a linear combination of 
the constraint set. Thus, for a primal posynomial program ' 
given by eqs. (3.25), (3.26) and (3.27) with all o’s positive, 
the corresponding surrogated program is 


minimize y^ (x) 


M 

subject to IZlYn ( 2 ) 

m=:l 

and x^^ > 0 } n = 1 , . , 


( 3 - 28 ) 

(3.29) 

(3.30) 


where nonnegative scalars that are normalized by 

requiring 



These will be referred to as surrogate multipliers . 
In order for the surrogated geometric program to be equivalent 
to the original program there must exist set of variables 
such that the optimal solution £ of Eqs. (3.28), (3.29) 
and ( 3 . 30 ) is optimal for the origina,! primal problem. For 
proof of the existence and uniqueness of the reader is 
referred to the book by Phillips and Beightler [58]. 


The solution procedure involves searching possible 
values of x £• At each iteration of the procedure 

a value of x will be chosen and the resultant S.G..P 
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problem solved. The solution that is obtained will then be 
checked for feasibility in the original problem and the 
chosen value for x will be altered such ths-t at each iteration 
it xvill better approximate the value of £. 

3.7.2 Problem Reframed in S.G-.P. --format ; 

The structure of the objective function remains the 
same while the constraints are reformulated as follox-irss 



Theroforo the one constraint in S.G.P. format can be 
represented as 






^i min 




74P+1 

^ax 




(3.32) 
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P 

Where (Yi + + Tzp+i + Y3PY1 + Yj-p+p) + Y4P.H = 1 

(3.33) 

P baso 3 Similarly the one constraint in S.&.P. format for 
Phase 3 can be represented as 




S. 

1 




y± s — ■ ■ " 

i=l i mas x=J 


S . . h f ^ . 

1 min , ■'v — ,, 

%+i S,. ■■ ‘ ^2H+i f„.. 


ai max 


H f . 

, ’c — ai min 


N 


^ai 


+ > Y/ 


Tl 


i_l 4b+i 


^6h+l — 
H 


ma: 


H ^ci H 

> (PZ **li ■'■ ill ®ik' i7 Yeppi J— 
i^ 1=1 ^ k=l 13. ■ n 


i=l 3 




“"ai 


B. . ^ri 


+ [ 


i=l 3=1 


Y . . ' “f‘ 

'^13 y 


fcl 


].<'! 


"'^ik V 


(3.34) 


wher'" 


B- . 
13 


/4- + t ) 

^\ii3 ci3 


3 

i 


!1L p i» * 

3- 0 w • 




%k 


^"'^mik "** '"cik 


haik 

“^ik 


) k = 1 




where 


h 

ZZ (7i + Y;,T+i + Y2]\T+i + '^3h+i ^.T+i 




i=l 


^ai 

+ P~ Y 


J^ri 

+ ZZ ik^ 761+1 


3=r ‘ k=i 


1 


(3.35) 
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Soluiiion Pr oce dure ; 

The steps that comprise the Surrogated G-eometric 
Programming Algorithm (SGPA) are as follows; 

1. Check if the problem can be formulated as a pos 3 niomlal 
and all 2 ^ are finite positive numbers. If not, SGPA 

cannot be applied. 

2. Determine if af least one constra.int is binding at 
optimality. If not the primal objective function can be 
treated as an unconstrained posynomial programming problem 
and solved. 

3 . Choose the initial values for the surrogate multipliers. 

In SGPA a sequence of posynomial problems having one 
constraint (the surrogate constraint) v/ill be solved. The 
coefficients of the terms in the constraint will be altered 
by cha.nging the value of x until an x" is located which 
satisfies the original constraints. The surrogate multi- 
pliers used in the problem that yields this are the 
required values of 

Obviously it would be beneficial to choose an initial 
value for X us close to as possible, thus minimizing the 
number of iterations reciuired by SGPA. If no information i s 
available with respect to which constraints are most likely 
to be binding, then it is reasonable to start with all 



Yjq = i/m, m = M. Using this starting va3-ue is 

simply an admission that there is presently no roasoii to 
assume any distinctions "between the constraints. 

If, however, some Icnowledge regs,rding the relative 
importance of the constraints is available in adva.nco, 
better starting Vfaluos for can be obtained. If a 
constraint, yj,^(x), is not binding, then = 0. Thus, if 
a constraint is not binding, it makes no contribution to 
defining the location of the optin^uii. If a constraint is 
thought to be particuls-rly important then its surrogate 
multiplier should be given a large value, while the surro- 
gate multiplier of a relatively insignificant constraint 
should be given a small initial value. 

4. Solve the surrogated problem using a suitable optimisation 
pro cedure . 

It should be noted that xj determined at this stage 
are the optimal values of x only for the particular program- 
ming problem obtained when a given value is assigned to 
X ,* is not necessarily the solution to the original 

posynomial program. In fact, these two quantities are equal 
only when x = . 

Each successive x xfritten as 

^ x^^* 'Ih.c if" 's corresponding to these x'® "'^'^ill 

be denoted by 
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5. Calculate the values of the constraint. 

Once the ve.lue of ^ has been determined for the 
(k+l)th choice of x? this value can be used to ascertain 
whether the constra,ints for the original problem 3,re sa,tis- 
fied for this If satisfies the original constraint 

set, then # a.nd = £" • If 2% does not satisfy 

the constraint set, then new values for x be selected. 

When 3^ does not sa.tisfy the constraint set, it will be 
convenient to subdivide this set into three mutually 
exclusive and totally exhaustive subsets, let 

0 = > 1 ] 

Q = (Si5;) -<1 ] 

= [r|yp (^iJ = 1 ] 

If ^ is the optimal value of the primal v.ariable 
vector for the origjjial geometric program, then when the 
constraint set is evaluated, 

0=6 (5.36) 

when this occurs the SGPA terminates. If Sq. (3.36) is not 
satisfied the algorithm continues to step 6. 

S. Alter the values of the surrogate multipliers.. 

If the value of lii corresponding to a x^ is such 
that 0 ^ 6f then x^ ^ a new value of x must be chosen. 

It is desirable to choose this new value ^ 



manner that yo*'~k+l^ — ). First calculate the following 
value , 


% = ^ (51 Yi ) y„ (2: )/ f: ra e 0 (3.37) 

Then calculate the now estimates of v usinv 


the follovjing equations 
k+1 k 


m 

kH"l 

m 

k+1 

m 


= U ''•n ■• 

= [H Yi - II Ol + 

IeQ ^ fep ^ 


m £ 0 
m £ R 


n ■ ^ (Yi - Yi'^^)]/n(Q)? 
i£R 


wnere 


A = a constant which nay be used to speed up the 


convorgence 

Y T = l^^-st v8-lue of surrogate nultiplior for which 

mXf 

7.-.(x) -^1, d £ R. ' 

iii 

Y r-r = last voJ-ue of surrogate nultiplior for 
'mil 

which yj,.^(x) >1, n £ R 

n(Q) = number of terms in the sot Q 

Elements of sot R a.re retained for 5 SGP iterations 
in order to avoid prenaturo shift to set Q. 

The multipliers are then normalized after equating 
to zero any multiplier with a nogativo value. 

' . , l.l . 1 fPlIR 
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7. Repeat Steps 4-6. 

Steps 4-6 are repeated until all the constraints 
are satisfied vrithin some predetermined tolerance. 

In the next chapter the results of the applica.tion 
of the solution methodology* given in this chapter, on s. 
test case are presented and discussed. 



CHAPTER IV 


RESULTS Am DISCUSSION 


The solution methodology suggested in Chapter III 
has been tested on a test problem available in literature [3]. 
The workpiece to be ma.chinGd on a multi-spindle automatic 
has the chara.ct eristics shown in Pig, 1, The various input 
parameters pertaining to the multi-spindle automatic and 
the workpiece a,re tabulated in Table 1. The results of the 
test problem are presented in Section 4.2. However, before 
we discuss the results, a. brief report regarding the computa- 
tional experiences and difficulties encountered in the 
implementation of the methodology proposed in Chapter III will 
presented. 

4.1 IMPLEMENTATIONS 

4.1.1 Difficulties Encountered ; 

Initially the Blau algorithm was directly used as a 
solution procedui-e for Phase 1 of the methodology. It was 
observed that the program was getting terminated due to 
excessive overflow conditions. The cause of this problem 
was not immediately traceable. It was decided to use the 
S.G.P. approach which reduced the Newton-ifephson matrix 
of the Blau algorithm from an order of 62 to 22, so that the 
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fciniG roauired for matrix inversions was grOf3,tly reduced. 

Moreover, the SGP approach rendered the problem no re tractable 

and eno.bled detection of the earlier source of error which 

was as follows s tho Blau algorithi'n roQuiros a 'good' initial 

solution to start the procedure, but since the initial 

solution used was not 'good' (refer Tablo 2), the procedure. 

for calculating the initial Lagrange multipliers, suggested 

by Blau, failed. This fe.ature caused trouble in tho S P 

approa.ch. too, at the initial stages. The source of error was 

detected after a, study of two test problems given by Dera.bo [13]? 

for which tho magnitudes of the Lagrange multipliers wore 

in the zangc of 0.2 to 2.0, while for the Phase 1 problem the 

magnitude was greater than 1000.0. Further, the conclusion 

that the procedufo for calculating the initial Lagrange 

multipliers was at fault, was reinforced i-rhen convorgencG was 

achieved after a.pproxinating tho Lagronage xnultiplior , for 

only the first iteration of S.G.P., by unity. This indicated 

n. X* x*3*ii*^ ©IB. 0'“'^ ij 

that by using this/tno first iteration of S.G.P. provided a 
•good' initial solution for the second, it oration. Though 
this arrangGEicnt provided convergence in Phase 3 also, it 
cannot be gonero-lised to extend to other problems. 

4.1.2 S. G.P. - a lote [ 58] ; 

It is intorosting to note that S.G.P. algoritln 
provides useful information at every itero,tion. This 
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information is not obta.ined whoii the problem is solved 
using other techniques. Notice that at the k-th it erections 
the problem, 

minimize y (z) 

0 

subject to 5 m = 1 ,..., M 

X > 0 

is solved. Forfinite k, it is the solution to a system that 
is the same as the original problem, given by eqs,(3.29)5 
( 3 . 50 ), except that the levels of resource availability 
have been changed. Hence it gives an optimal solution to a 
model that might be an acceptable alternative to the 
original model. 

Another important by product of the SGP algorithm is 
that it provides information for sensitivity analysis. At 
any iteration the vali;es of the surrogate multipliers for 
the problem that has been solved at that iteration are known. 
In particular, when the optimal value of the original 
program is found, .the optimal values of the surrogate multi- 
pliers are known. 'These multipliers give the relative amount 
by which the optimal value of the objective function would 
change due to change in resource levels. Hence, the 
surrogate multipliers indicate which constraints are most 
binding and which constraints are of least importairce. 
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Lot UG now turn our attention to Table 3. Prosonted 
in this table is the progress of the optimizo.tion algorithm 
(SGP) in terms of tho constraints (which w’ore tight or nea-rly 
tight at termination of tho progra-n), the corresponding 
surrogo-te multipliers and the objective function, for 
Pha,Go 1. It has. been presented to indicate the efficacy of 
the updating procedure for tho surroga,te multipliers, which 
has bot.;n developed for this thesis (a procedure suggested 
by Stao.ts in [58] was found to proceed at a crawling pace). 
The first column of this table indicates the number of 
iterations of the Blau a.lgorithr."'. required for the corres- 


ponding S.G.P. iteration. The nomenclature for the symbols 
used is as follows; 

yQ is the magnitude of tho objection function 

y^,..., y2_Q are the magnitudes of the upper feed limit : 

constraints for oporations 1 to 10 rospoctively i 
^11’**’ ^20 magnitudes of the lower feed limit | 

constraints for oporations 1 to 10 resp. i 

y^i,.., y^Q are tho magnitudes of tho upper speed limit 

constraints for oporations 1 to 10 resp. | 

are tho nagnitudos of tho lower spood limit j 
constraints for oporations 1 to 10 rosp. J 

is the magnitude of tho power constraint 


^ 31 ’ **> ^40 


y 


41 


and Y- is tho surrogato multiplier for constraint 


y 


1 »■ 


1,2, ... , 41 > 
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Further e. constraint value 

= 1 ind rentes the.t the constire-int is tight 

•<;; 1 indicates that the constraint is loose 

>1 indicates that the constraint is \dLol.ated 

Based on the 8 ihove discussion on S.G-,P., referring to 
Talole 3 we find that at the 8 th iteration of the S.G.P. 
lalgorithra, the power constraint is the most hinding one 
(T 4 I = 0.835 while the values for the other nultipli-rs o-ro 
below 0 , 02 ). Further the value of the surrogate multiplier 
indicates tho.t the sensitivity of the objective function tc 
the power constraint rola.tivo to the remaining constraints 
is 88.5 percent. 

Now we shrrll discuss the results of the t est problem. 

4.2 EESULTS,; 

4.2.1 Pliase^l; 

The test problem has 10 operations , This respl.ts in 
20 variables and 41 constraints. Corresponding to each 
operation there would be one variable j;or speed o.nd aeiother 
one for food. For each variable there are two constraints 
corresponding to its upper raid lower limits. Further, there 
would bo a constraint resulting from the total power available 


on the machine tool. 



laolo 2 presents tlie initial solution solectod for 
fch.o cutting po.raEiGt ers , The initial food valuos for Go.ch. 
oporation WGro choson within the limits provided by tlio 
Taylorized rr.nge. liowovor, the initial valuos for speeds 
tiTGrG cliosGn noo.r tho lower spood limits, decided by the 


Ta.yloriZGd rc0.ngG. This wa.s done to avoid the violation of 
the power constro-int. Table 2 also presents the optiraa.1 ssj 
VcOlues of the cutting parameters obtained, for ea-ch oporcation. 
It should be noted, with rofei'ence to Table 3> tlmt for tho 
optimal solution obtained two constraints are getting violo.ted 
thos.,:; constraints are; 

1. the fe.d constra.int for opero.tion S. As y^ is 1.003, 

the violation of the constro.int is hy 0.5 percent. 


2. the power constraint. As y^^. 1.05, which imp3.ics 

that tho violation of the constro.int is b;=" 5 percent. 

Those violo.tions were not considered excessive. As a 
matter of fo.ct, the criterion shop tod T'Jp.s that s.ll tho 
constraint violation should be loss than 5 porcent. 

The Phase 1 took CPU time of 4 minutes on IBM 7044. 

4.2.2 Php^e,,2y 

'.Uho output of Phase 1 is used as the input for this 
phase. Tho feature vectors of the varioois opero.tions are 
prcsontvcd in Table 4. The results of the clustering are 


given in 


Tables 5(P) 5(0). 



fable 5(B) gives the results of grouxoing ulicii only 
ch .0 j..iiriii.iui.i and ittaxiiiiUiU. speed o.nd foed lir-its I’.nd ■Di'!;CGdonce 
vialucs ,are considered in tlio fooefcuro voctori in tlie computer 
progran this was acconiplishod by equating a,ll the optimal 
speed and feed vsilues of the foaturo vectors to zero before 
grouping. Tlio inultixDlication factor for the prscedoiice va,lues , 
requirod to just satisf^r the p3x;cGdencc constraints, is 1.22. 

Table 5(C) gives the results of grouping when onljr tho 
optimfxl feeds and speeds (of Phase 1) and precodonca values 
n,ro consideredj this vras accomplished in tlio computer program 
bx?- eq_uating the remaining elements of tho foa-turc vectors 
to zero before grouping. The required value of multiplication 
fo.ctor is 0.74. 

The CPU tine required on IHIi 7044 in each case was 
0.5 minutes. 

4.2.3 Pha s_e , 34 


Tho grouping; information obtained in Phase 2 considorin 
optimal speeds and feeds and procodence values was used as 
input for this phase. Tho problem was formulated for a 
3-sp indie automatic which yielded S va.riables - 3 speeds, 

3 a,xi;,il foods, sind 3 rs,dio:.l feeds — and 29 cons tn-ixicu 
1 power constro.int, 10 cycle tino constro-ints , ano lo upper 
a.nd lower limits for the 9 rr-ri-ables . Tho conpubo.. 
progro.m was set to toiaiinate when violation of all constreinti 





60 


is loelov^ 5 
paraniotei’S 
To.blc 6. 


porcGiio, TiiC’ rcsul'ts, op'cinial v.aluc8 of 
and tlio initial solution, arc presontGcl 
Tlio initioJ. values of the cutting po-ranct 


cutt ing 
in 


for a. group wore chosen o.rhitrarily frora aiaong the Phase 1 
optimal values of tha/fc group. 'Tlio C 3 rcle time ohtainecl was 

0. 93. rainutos. The constraints violated and tho extent of 
violation aro listed helows 

1 . Upper food linit constraints for operations 3 
and 9 woro vio 3 -atod hy approxiiTia,tely 0.3 percent, 
raid that fox' operation 2 hj 1 percent. But after 
rounding off the corresponding feed vario.blcs to 
2 decimal places the constraints bocano tight. 

2. The pov'r constraint was violated by 4 percent. 

The CPU tine, on IBM 7044, requirod for this phase 

was 2.5 ninutos. 

4 .2.4 Ccmpari.§!il-l„i?.^. 


'Two sots of results obtained by Barto.lucci ct -\1 [3 
port.i,ining to grouping of operations on'o pres one ad in 


Table 5(A) 
spindlos ) , 
occur tog 0 ' 


. Those results perto-in to 3 groups (or 
It is obsorved that the following opere.ti 
thor in both sots, 


1. Operations 3, 7 ?.nd 10 

2. Oporr.tions 2, 4 and 6 

3 . Operations 8 and 9 

4 . Opomtions 1 and 5 
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Taking these 4 groups of operations as the basis for 
comparison, we conclude that case C, i.e., the case in which 
only txie opbimal speeds and feeds (of phase 1) along with 
precedence relationships are considered, approxima.tes the 
results of Barta^lucci et al mors closely. 

The opti'nization results of phase 3 give a cycle time 
of 0.93 minutes and the total CPU time, on IBM 7044, req.uired 
for all the three phases is 7 minutes. The corresponding 
results obtained by Bartalucci et al are presented in Table 7. 
Prom this table we observe that the cycle time of 0.93 minutes 
obtained in this study is ver 3 '' much less than the results ob- 
tained by Bartalucci et al. Moreover, the time required for 
the computation, using the 3-phase procedure, is much lower. 
Further, it should be noted that IBM 7090 (used by Bartaluccd 
et al) is computationally faster than IBM 7044 by a fa.ctor of 
1,4. 

Table 8 represents the allocation of groups to the 
spindles of the 3-spindle automa.tic along^-irith the type of 
tool and slide type for each operation. 

4.3 SCOPE FOR FUTURE WOEK; 

The present work considered the determination of 
optimal cutting parameters for multi-spindle automatics 
without giving due emphasis on the evaluation and selection 
of alternate tool— change policies. However, it is a well- 
known fact that optimal cutting parameters and the choice 
of todl replacement policy are related, therefore, it 

would be worthwhile to develop models for the determination 
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of optiEir.1 tool roplacenont policy md optinr,! cutting 
paramotors in oji integrated fashion. 

In th-G presGiit work we liavo assunod tho classical 5 
d eov^'ri-i inis tic Taylor's tool life equation. Howovor, raany 
investigators Imve roportod in literature that tool life 
is a probabilistic phenonenon. Various typos of distribution 
3_iko ''Jcibulj, log-noma.! ha,vc been reported in litei’-ature 
[ 46 , 51 , 52 , 53 , 62, 63 ] for various tjrpes of tools. Vurtlicra 


I'riodnan and Tipnis [25] and Ravigiani et al [55] have 
prosoiitod a. nox^ perspective, vis,, the Cutting-ratvi - Tool-life 
(R-T) characteristic curve x^jhich represents tho loci of all 
optimal combinations of speed, feed and depth of cut, i.e., 
cutting rate, and tho corresponding tool life valuos. The 
potential of this curve with respect to cutting parameter 


optiraiization norits investigation. 


A new trend in manufacturing called Group Technology 
is assuming the proportions of a break-.through and opens up 
a whole range of optimization problems for nulti-spindle 
automatics. Group Technology could be introduced as a technique 


for identifying and bringing together related 
components in a production process in order to 


or similar 
take advantage 


of their similarities by making use of, for exampl-C, th<= 
inherent economies of flow production methods. ihus one needs 
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to develop ''lodels and Rothodologios for the adaptation 
of laulti-sp India a.utono,tics to "batch production environ 
nonts considering Group Tochnologjr concopts. 
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Ttuble 2; Input /Output of OptimizntionsPhase 1. 




Initial guess 

Optiinun "W-lues 

Process 

Typv; 

Opr. 

Ro. 

Speed 

(rpm) 

Peed 

(imn/r:’V. ) 

Speed 
( rpn ) 

Peed 
(mm/ rev. 

— ing 
. time 
) (Optimo.: 

Radial 

1. 

60.0 

0.25 

lOS.O 

0.4 

0.59 

Axial 

, 2 

60.0 

0.2 

123.0 

0.37 

0.99 

Radial 

3 

60.0 

0.1 

215.0 

0.09 

0.78 

Axial 

4 

SO.O 

0.35 

70.6 

0.58 

0.98 

Radial 

5 

60.0 

0. 25 

85 . 6 

0.4 

0.73 

Axic'.l 

6 

70.0 

0.35 

111.5 

0.54 

0.76 

Radio.l 

7 

60.0 

0.1 

99.0 

0.19 

0.79 

Rad in.l 

8 

120.0 

0.25 

198.6 

0.35 

0.5 

Axial 

9 

120.0 

0.25 

235.0 

0.5 

0.8 

Axial 

10 

120.0 

0.2 

237.0 

0.18 

0.23 

All va 

lUGS 

of vspeeds 

are round ed 

off to 

one deci 

,mal place 

All va 
two de 

lues 

cinial 

of feeds ca 
places. 

nd pro cess in 

g ti®es 

are rounded off t 



Table 3: Progress of Optimization Algorithm (SGP); Phase 
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Table 5° Operation Groups. 

A; Groups obtained by Bartalucci et al [3] 

Group No. 1 2 3 

3 Groups with 3,7,10 2,4,6 1,5, 3, 9 

40 HP 

3 Groups with 3,7,10 2, 4, 6, 1,5 8,9 

60 HP 


B; Groups obtained using only maxima, minima of speeds and 
foods and procodence values of feature vector. Hulti- 
plication factor neccsson-y for precedence values = 1.22. 


3 Groups with 3,7 2,6.9, 8,4,10 

40 HP 1,5 

C; Groups obtained using only optimal speeds and feeds and 
procodence values of feature vector. Multiplication 
factor necessary for precedonce va].ues = 0.74. 

3 Groups with 3,7,10 8, 2,6,9 1,5,4 

' 40 HP 


Values indicated in the table e^re operation 


numbers labelled as in Pig. 1. 
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Table 7; 

Piosults obtained 

by Bo.rtalucci et nl 

Humber 

of 

spindles 

Power 

(HP) 

Cycle time 
(min. ) 

Comput-ation time 
on IBM 7090 
(min. ) 

5 

60 

1.51 

28 

3 

60 

1.54 

8 

3 

40 

2.42 

8 
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Table 8; Allocations of Groups to Spindles. 


Station Group 

No, Allocated 

Operations 

Tool type 

Slid e type 


2 

Overbanging 
turning tool 

Axial 

1 2 

6 

Turning tool 

Axial 


9 

Turning tool 

Axial 


8 

lacing tool 

Radial 


10 

Recessing tool 

Axial 

2 i 

3 

Flat form tool 

kdial 


7 

Flat form tool 

Radial 


4 

Turning tool 

Axial 

3 3 

1 

Facing tool 

Radial 


5 

Facing tool 

ifeidial 


dll 



All dimensions are in millimeters 


Pig.1 : Test Problem ; Job to be machined 
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APPEIDIX A 


BLAU (CtOMTRY) AI&ORITHi 


This algoritlani finds the nininun of a nultivariahlc , 
nonlinear function of geometric form represented Td],?- eqns. 

( 5 . 25 ) V ( 3 . 26 ) and (3.27). The particular algorithm included 
in this section was developed by Blau [5] and programmed 
by Garcia and Hogg [31]. Thu notations used are t ho same as 
for oqns. (3.25) to (3.27). The method proceeds as follows; 

1 ) Enter problem with input parameters. 

2) Determine initial weights; 


T 

t=i 


ot '^ot 


TI 


'oxn 


n=l 


n 


Dual function, V 


"nit 


Z 


K 

TT 


a, 


mtn 


n 


m = 0,1,. 


? -i- y • • • S? 


M 


3) 


■" n=l 

^ot “ ^ot^'^ 

Calculate the vector of orthogonal conditions for the 


objcctivo, and the matrix of those conditions for the cons- 
traints : 


T. 


m 



A-2 


~ ^ot n - I5 . . . I'l 

^ ^ .yvaD-uatc; initial Lagrange multipliers 

X = K)~^ ic'^ H (IC*' denotes K transpose) 

5) If this is the first iteration, go to Step 6),* 

o thv.irwis 0 , dotormino now ¥eights as in Step 2), neu orthogonal 
conditions as in Stop 3), and new multipliers, as follows j 


^no 


w 


“Old 


f>) 


low proceed to Stop 6). 
Calculate the matrix T 


tj ^ 


T 
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[ m ‘^mt ^mti ^mtj l^mt^ 


T. 


t=l 

Evaluate error 
In 


X--- ^ot ®'mti ^mtj ^ot 


YZ ®ot =‘otl ^ot 
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M m 


jL — _L 5 * • » 5 
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- i: [ i: “mt V 
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o Ot 
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T. 
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'^mt ^mt 


1 , . . . , N 


i = LT + 1 


N+l+m ; 
m=T, . . ,M 


1 



A^-3 


1 2 


N N+1 li+2 


H+M+1 


1 1 

2 2 

n: 

N i] 

J N-hl 0 


11 
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N+2 1 

N+3 m 2 

• • 

# • 

N+Mh-1 M 
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■T 



9) Invert the matrix R. 

10) Rind vector of adjustments; 


11 ) 


R 


iiln X 
Ain V 
Ain y 


Calculate new values of independent variables: 


S = z exp (Ain x) 
V = V exp (Ain Y ) 


12) Has solution converged to acceptable limit? 

Yos } print results and stop. 

No; go to Step 13. 

13) Has maximum allowable iterations been reached? 

Ygs_j stop, print message. • 

Hoi return to Step 5 • 

A flow diagram illustrating the above procedure is shown in 


next page. 















APPENDIX B 


SURROGATE GEOMETPJC PROGRAM 


Usage ■; 

The program consists of a, main program and throG 
subroutines, GP3, GPIO, and GP22. The main program roads 
the data, prints the results, and performs a rna.jor pa,rt of 
the algorithm. 

Subr outin es Roqu irod ; 

Subroutine GPp calculates the absolute value- of all 
terms in the problem. Subroutine GPIO inverts the Newton- 
Rap hs on matrix, or ono of its submat rices. Subroutine GP22 
calculates orthogonal conditions for the objective function 
and the constraints. All linlcago between main program and 
the throo subroutines is through COMMON statements. 

Dos e rint ion of Parameters ; 


N 

Total number of variables 


M 

Number of constraints 


NO 

Number of constraints in the original 

Geometric 


Program 


TS 

Assumed sign of the optimal value of 

the objective 


'function 


MRX 

Dimension of raa trix Rj (M+N+1 ) 


CONVG 

Convergence criterion 




3-2 


J\i Vector giving nunibor of torras per polynoiaial 

G Voctor of coefficients 

A Array of exponents 

A Voctor of initial estitnates of optimal solution 

¥ Vector of weights (dual variables) 

E Error voctor 

AMD Vector of Lagrange multipliers 

PIVOT Voctor used in matrix inversion 

LY Card reader unit number 

LZ Printer unit -number 

SIM Objective function 

V Dual function 

KTT Total number of terms in the problom 

Ml Total number of polynomials in the original 

Goometric Program. 

SUE Voctor of surrogate multipliers 

DIMENS lO I Roquiromon ts ; 

The DIMENSION statements should be modified according 
to tlio roquiromonts of each particular problem. The para- 
meters in parentheses included in the following statements 
conform to the Parameter definitions above; 

DIMENSION CCETT), X(N), ET(M1), ¥(KTT), A(J;2rT,N)s 

R(MEX; MEX), PIVOT (MEX) , AMD (M), E(S:TT), 
IPIVOT- (MEJC), INDEX (MEX.„ 2), GO (KTT), 
SUE (NO) 



B-3 


DTJIEKS ION 


Imout Formats : 
Card typo 


NC), WEI (EC), Kin (EG), M(EC), 
FlA(jIO), lA(NC), SO(FC), STO(l'TC) 


EOK'IAT 


COEEEHTS 


1 (214, E5.1) I". M, TS 

2 (214) • NC, MPJ[ 

3 (14) liT(L), L = 1, Ml 

(One typo 3 card for oacli polynomial, including 

objcctivo function) 

(5E 15.3) C(I) (Coefficient Itli terra) 

(If I is greater than 5, additional Type 4 cards arc 
roquirod) 

5 (5B 15.3) A(l,j) (Exponent for J-a^ 

varisibl^ I'bli ‘tGrin) 

(If J IS sroator tten 5, aldltional Typo 5 cards roquirod 
for each term.) 

6 (P 11.5) x(J), J = 1, tr 

(Ono t;-p= 6 card requlrod for each rarlallo) 

Out pu t J. 

Tho problem statistics which include, numbsi 
variables, number of oonstramits. number of terms, and 
aign Of tbo obioctiuo function are printed. Tberoaftor.on 
each S.S.P. iteration, tbo current solution, absolute 


3C 
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multipliers arc printod. This is continued until the 

convergence criterion hc?.s hoon' satisfied, after which tho 

message 'CONVERGENCE CRITERION SATISEIED' is printed. 

Summary of User Re q ui r ements j 

1. Formulate problem according to requirements 
explained in Section 3.7. 

2. Specify values for N, M, TS, NC , MRiL, LZ, lx, COlWQ. 

3. Specify initial solution, X(J). 

4. Adjust DIMENSION statements in main program and 
subroutines. 

5. Modify FORMAT statements 1000, 1070, 1100 and their 
associated R:f;AD and WRITE statements for the particular 
problem being solved. 



APPENDII C 


PROGRA.M PO.R GROUPING OP OPERATIONS 


Usa^.; 

Tho pi'ogram consists of a main program and five 
subroutines BUILD, lUffiAIS, SINGLE, MIN and FiAZ. The main 
program reads the data and prints the results. 

Subroutines Requir ed : 

Subroutine MIN finds the minimum of a set of values. 
Subroutine MX finds the maximum of a set of values. 
Subroutine SINGLE updates the statistics for each -'group. 
Subroutine IfflEANS determines the closest group for each 
operation. Subroutine BUILD initializes the statistics of 
each group and calls KlIEANS. All linlcage between the main 
program and the subroutines is through ar unients of the CALL 
statements. 

De scription of Par ai' 3.e t e r s ; 

K Number of spindles 

M Number of operations 

N Number of characteristics for each operation 

ITER Number of iterations for convergence of the Mac Queen 

procedure 

ENC Increment or decrement for the divisor which 
adjusts the precedence values (-1 ENC < l) 



Vector of mrcciraura values of all characteristics 
of niojCLi-iuiii values of a3_l characteristics 
for radial oido rations 

VJA Vector of raa,simum values of a.11 characteristics 

for axial opere.tions. 

A Array of chaiact eristics for all opera, tions 

TIM Vector of cutting tiraes for all operations 

OLL> Vector of cutting lengths of all operations 

i;’A.S Vector of consisting of binary values » 0 indicates 
a radial operation and 1 indicates an axial 
operation 

I'TAS Vector of operation nunabors (corresponding to Fig-l) 
DOLUS Vector of distances from group centroid 
DIMEES IO.I Requi rement s; 

The DIMElTSIOl statements should be modified according 
to tho roquirements of each particular problem. The parame- 
ters in parentheses included in the following statements 
conform to the parameter definitions above; 

DIMEUSIOI WR(N), UA(U), TIH(M) , 01T(M)j FD(E5M)5 
MClU (K,H), ¥(I), SUM (8,U,K), 

3[(1T), ITCLUS(M), DOLUS (M) 

M(M), MS(M), M(M), HAS(M), DX(M), 

JGLUS (M), X1(H), X2(M), X3(M), X4(M) ,lTST(ii) 


DIMSUSIOI 



Inpu t F orraats: 


Card 

Type 


Format 


1 (7 F12.4A3X 

I 1, 13 ) 


Contents 


(A(I,J), J=1,]\T), TII'Ul), OLT(I) ,FiAS(l) , 
MS (I) 

(Values corresponding to tlie J'tF _ 
cliarac tori Stic of tlie xth opercvfcion) 

(One card is recgaired for the Itli 
operation), I = 1, H) 

w(j), J = 1, F 
im(J), J = 1, 11 
VA(J), J = 1, 11 

EITC 


Output ; 

Tho information in the print-out includes tho ^-roao 
numbers, members of tlic group and maximum and minimum speed 
and feed limits for each group. The raultiplication factor 
for tho precedence values is also printed. All this infor- 
mation is repeated for 5 different values of multiplication 

factor. 

S ummary of 

1. Formulate pro b lore according to requirements explained 

'in Section 3.3.3. 

2. specify vnluop for K, M, N, I'lBE In tho BATA 

st at Gill out . 



Mult ipli f ic at i o n facto r is given by [ 1 / ( ' ■ ( 6 } i ) J 
so specify val^^e of EEC so as to obtain roquirou 


range for pirocodence valu s. 

Ad oust DIMEESION statements in ma-in program and 
sub rout in es , 


Modify FORiM-T statement 2 and the associaued BEAD 
statements for the particular problem being 


solved. 
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!:*««*« NUMRKP OF TFF;1«S IM EACH FOtYNOH'IAL 
KTT s 7 v. 

PEAnClY, UlBtrO CKTcr), X«J,Mi5 
#*«## iniruillt eUPPOGATF CONSTANTS 
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SIG 
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15 
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r« X j» »* 1 '''r.''-'’^- '> TPK ,'„rfjKL 

f l' M C ('n'jl,Xf2in,E1C'^),V.(7yn,AC7t;,2c-'),Pf«7,22)»HVLTC7?), 

1 ^ 'U M .E f 7;o 

r O-'i- r'f V ,r J VU7 ,E,P.SUF,V,SIC:,AE'f:,C.XfKT,A,i*f’ ,?'• .•''f » ?■ TQ.f i- , 
t ■ 1 r !T ,»''b ! ,r-K 

' ' ' 1 '7 5 \ 

1 >1 17^' t* 1 f ^'1 

M s 12 + 1 ' ' ^ : '7 ' 

? V = 1 1 4 f 7 ( L ) - 1 

I'O 'pi I = U , T2 - V . ■ ■ 

CM s 1. . ^ 7 ^ ^ ^ 

f'O TEv* .7; ^: ' ; ; 7 

HM 5 XUD 

n (Hful ,E‘n, 0,e) .GO TO., 462 - . 7 „ ■ : 7 . -^ ' 

r 1 s A ( J. » 0 ? 

fAT ’,F-(/, GO 10 46 vT 

r-BJ a FD0*<»AT ., . . ... 

> ( .f ) s it.(l) * F-BJ 
4 60 ri^fTTBUF 

OCi > ( n s « AP3(C(I)) 

FETUPN 

E ND ^ ^ - ■ ' ■ 

C*»^!^^^*^t^»«t#»«»l#•■t^»«■»*^‘'»«**•»«'‘^■******■•*■*^'*''**'*‘*^**''‘^******'*‘'^*^**'****''’‘^**'*'** 



r»»i»««'THF SIIPBCUT IDE GP2? IS USED Tq FTNP THE VECTOF CE OPThOGCDAl 
C#«<HH*rODDlTIODS FOP. THE OBJECTIVE* ADC THE MATPIX OF THESE CGKCITIONS 
C#4Ht«»tF OP ’ the C0H5TPA1DTS 

DIHEDSIOP C(705,XC, 20)*!<TC2)f'''«C72)fAC70,223*fi(22*22)»PIVCTC221i 

‘cO«MO^'^i’?iV05,E,P,5Ut-,v,5IG,OM,C,X,l<l,*,l«K.N,«',l'l.l‘TT,KTO,«B, 

iWl# IT*KB} ,HP 

DO 4e@ iPslfT^ 

j 2 « 0 ■ ■ ■■■ '■ ' - ' ■ ' 

DO 4J0 JP»MI,MP 

Ij » OH • M 4 1 

II « 12 4 1 

T2 » It 4 FTCL3 »• 1 

a « IB . 

DO 400 

480 PClR*jP)'^*^HtI^rJB)’'' 4 CCXI ^ A8S(C£I)) « PJJ * F ^ 

PEtUPN • ^ '■■ 

END' 


APPENDIX C 


PROGRAM FOR GROUPING OP OPERATIONS 


Usa,'5:e 

Tlio prograrii consists of a main program and five 
sulD routines BUILD, JUDEANS, SINGLE, MIN and MAX. The main 
program reads the data and prints the results. 

Sub rout in os Requir ed : 

Subroutine MIN finds the minimum of a set of values. 
Subroutine MX finds the maximum of a set of values. 
Subroutine SINGLE updates the statistics for each group. 
Subroutine ICMEANS detenTiines the. closest group for each 
operation. Subroutine BUILD initializes the statistics of 
each group and calls MEANS. All linlcage between t he main 
program and the subroutines is through ar uments of the GALL 
statemonts. 

.DGOcript io n of Paraiaeters ; . 

K Number of spindles 

M Number of operations 

N Number of characteristics for each operation 

ITER Number of iterat iongfor convergence of the Mac Queen 

procedure 

ENG Increment or decrement for the divisor which 
adjusts the procodence values (-1 < ENG < 1) 



'■f Tcctor of nmsimum values of all charactorist ics 

!■ .R Voctor of rua^cLniuin values of all characferistics 

for radial operations 

"doctor of raasimum values of all characteristics 
for axial opera-tions. 

A Array of chaiact eristics for all operations 

TIM Vector of cutting times for all operations 

OLD Vector of cutting lengths of all operations 

I IAS Vector of consisting of binary values « 0 indicates 
a rfxdial opera-tion and 1 indicates an axial 
operation ■ 

NAS Vector of operation numbers (corresponding to Fig.l) 
DCLTJS Vector of distances from group centroid 
DIMENSION R eq u i re m ent s : 

The DIMENSION statements should be modified according 
to tho requirements of each particular problem. The parame- 
tora in parentheses included in the f oil ov^ing statements 
conform to tho parameter definitions above; 

DIMENSION WR(N), NA(N), TIII(M), OLT(M), PD(K,M)j 
MOLD (E,M), A(H,N), ¥(N); SUIi (8,N,lO, 

X(N), NCLUS(M), DOLUS (M) 

min), MAS(M), NA(M), NAS(M), B{M,N), DX(M), 
JCLUS (M) , X1(M) , X2 (M) , X5 (M) , X4 (M) , N3T (M) 


DIMENSION 



Input Porraats; 


Card 

Type 


Format 


(7 P12. 4,132 
I 1, 15 ) 


Contents 

(A(I,J), TIH(I), OLT(I),mS(l) 

I'TAS ( I ) 

(Yaluos corresponding to the Jth 
cliaractoristic of the Ith opere/fcion) 

(One card is reo,uired for the Ith 
operation j I = 1, M) 

V/( J) , J = 1, I'l 
WR( J) , J = 1, IT 
WA(J), J = 1, H 


The information in the print-out includes the group 
numhers, members of tho group and maximum and minimurii speed 
and feed limits for each group. The multiplication factor 
for tho procedence values is also printod. All this infor- 
mation is repeated for 5 different values of multiplication 
factor. 

Su mmary of User RequiromentrS ; 

1. Formulate problem according to requirements explained 
in Section 5 . 3 . 5 * 

2. Specify values for K, M, 2, ITSR in the DATA 


statement . 



C-4 


Multiplification factor is given by [l/(^'(6) - f 
so spec if 5’’ value of SNC so as to obtain required 


range for precedence valu 'S. 

Adjust DII'IBNSION statements in main program and 
subroutines. 


Modifjr PORMT statement 2 and the associaced EnAD 
statements for the particular problem being 


solved. 



C#«»#I* rpOGPA?' rOR rtjf'STKPlKG OPSPATIONS ONTO THE SPINDI/PS AVfllPAHI.F 
nTMFNSION WR(73,lAiAC7),T.tMClP],Or,T(Jf?),HU&it^?)»WclUf5,ie3 

nTMFNSICN A(lF, 73 ,W( 7 ),SUP(R, 7 , 53 ,X( 75 ,NCTlJStie 3 ,DCH;SCie 3 
nTMPPSXON MA(1(41,FASflB),FACie!),NA5(l(n,P(lt?,73,D)f(Jg) 

UTMR^SlC^ OritlS(10),Xl(l^l,X2f Un»X3(1tJ)*X4ClP.3»NST(lF3 
1 FTR^'AT (7F12,4,2I9) 

7 FORMAT (9F7,3,UX,I1,I33 

3 FORFAT ( 7 X,I 4 , 9 H GLUSTEFS ) 

4 FrR^•A7 ( 10X, HHftUSTIF RliFBEB ,14) 

FOH^'AT (21X,1?19) 

6 FOR^'AT (/21X,UHFEED RADIAL , 4HHINs: , F7 , 4 , 4 X, 4H« AX= ,F7 , 4/ 
t 2lX,l0FFEEn AXIAL ,IX, 4HKINa ,F7.4,4X,4HFAXs , F7,4/ 

7 2\t, 5HSFEED ,6X, 4HMINs ,F7 , 1 , 4X , 4HMAX= ,F7,l//3 

7 FnRMATC21X,T9,FU,UFJ3.2#Fie,3) 

« FCFMATC27X,4HCPFN,UH SPEED ( PPM 3 , 1 3H FEED CMM/REV 3 , TIME (SEC 3, 5H 

C H,P,/3 

7 FCPMATC21X,39W MULTIPLICATION FACTOR FOP PRECEDENCE a ,F6,2/3 
DATA K/5/,M/tC/,N/7/,ITEB/3/ 

READ C5,2)(CAfI,J3»a=J |N),TIM(I3,OLTCI3rMAS(I), NASCI5i I = lfM3 
PFAD (5,23 (WCI3, 1*1, N) 

PFAD C5,23CWA(I3,Tsl,N) 

RFAD (5,23(WB(n,lel,N3 
PFAD (5,2) FNC 

C*#*HK* STORE IN ORDER OF RADIAL FIRST AND IREN AXIAL CPERATICN^^^^^ 

T, B 0 

MO at 0 ; 


no 15 Ujsl,2 

IF (Jil ,EQ. 2 3 MO » 1 

LM * L 

no 14 1*1, M 

IF (MAS(I3 .NF, M03 GO TO 14 
L * L + I 
DO 13 Jsl,N 
HCL,J) « A(T,J3 

13 CONTINUE 
HA(L3 « MAS(13 
NA(L3 » NAS(13 

14 CONTINUE 

15 CONTINUE 

STORE PRECEDENCE VECTOR ■ 


DO 20 l*l|M 

ce^n'^aBa,.,, 0 * 1 , 7 , . «cn, nmi,, i-i.ib, 

?::::: iziiiii «th 

CALL MAXCTIM,M,TIMAX) 

DO 22 J»lfT 

D(3 22 I«lfLM 

22 nT«u JYTST MIYTMA 

C«#i#*N 0 RMALl 5 E AXIAL OPERATIONS WITH AXIAL MAXIMA 


LMl*LM 4 rl 

DC 23 a«l #7 



WvIsWAfU!! 

DC 2 3 

23 P(T»a)=p(T,.n/w\j 
nc 25 1 = 1 , M 
m sctTci) 

BT7sBCI,7) 

SFEFP = 0I,/(TTMAX#B17^ 

IF (RPEEr.nT.RfI,3)3 SPEED=BfI.3) 

IF(SPKfc;r.l.T,BCi#4)3 SPEFD=E(I.4) 

TIMCI)sOL/CPI7#SPEED) 

BfI,5)sSPEED 
25 COMTIMUE 

DO 200 MJsl ,5 

CHANGE WEIGHTS rCP PRECKPENOr 
DC 28 T = l,1,0' 

28 BCI,6) = DXCI1/W(63 
W(6) = WC63 - ENC 
no 200 KK=1,K 

If (EK.EQ.l .OP, KK.EG,43 GO TO 200 
MB s 1 
M s LH 

N = 7 . ■ . . 

NN a 1 
K3 a 0 ^ 

THE NEXT FOUR PINES APE USED WHEN ONLY OPTIMA APE CONSIDERED 
no 2<3 1 = 1,10 
DC 29 Jal,4 
ACI,a)8RtI,a) 

29 

30 CALL BlJlLnCB,M,N,K,ITEP,SUM,NCLUS,DCL0S,X,NA,NN,MM,KK3 
no 35 1=5,7 

DO 35 K6=t,KK 
KH6 a + K3 
AfKK6,I) a SUM(3,I,K6) 

35 CONTINUE 

IF (M ,GT, LM3 GO TO 40 
M a 10 ■ 

K3 * FK 
MM a LM + 1 
GO TO 30 ■ 

40 CONTINUE 

THE NEXT THREE LINES ABE USED WHEN ONLY OPTIMA ARE CONSIDERED 
DC 50 Ial,10 
DO 50 Jal,4 
50 BCI,a)aA(I,J) 

C«#«#« STORE IN A matrix THE NAX/MIN VALUES OF SEED AND FEED FOR 

THE RADIAL AND AXIAL GROUPS TO PREPARE FOP PADIAL/AXIAL GPOUPIFj 
K1 0 1 ' 

K2 R LM 
K 3 = 0 

60 DO 80 KMb1,KK 

K4 » KM + K3 . . 

L a 0 ^ . 

DO 70 I«Kl,K2 

IF (NCLUSd) ,NE,. FM3 GO TO 10 



NSTCn S K4 

1. = 1 + 1 

XI fT.l a B(I,n 

X.1(L) a B(l , 3) 

Xl{!.3 s 
7«* cnNTIMlI^ 

fAl.L MTNfXl,l,ACK4,U3 
tat.!.. MlN{X3»l#A(K4,3)) 

TAIL MAXCX2,I ,A(K4,2)) 

CALL MAX(X4,L,AtK4,4)) 

80 CONTINUF 

IF (K3 ,EQ, KK) GO 10 90 
K3 a K'K 
K1 s + I 
F 2 s 1 0 
GO 10 90 
90 CON'’'P'UF 

PEOCPMALTSE with RFSPECT to OVFRALt ^aXIHA 
KK2*KK# 2 
nc 95 J»l»4 
WPjsWPtO) 

WJstW(JJ 

nc 95 I*lrKK2 : 

IF(T,GT,KK) WPOsWACJ) 

95 A(I»J)»A(T,v1)#WPOI/wa 
NN a 3 


Ns 6 

105 CONTINUE 

M ® K 4 . 

y w ssa 1 

DC 100 l«lpK4 


100 StJh RADTAL/AXIAn GPOUPING ON BASIS OF SPEED AND PPECEDENCj 

^rNExi UorUNEf IpE used IhEN only optima ABE CCNSIDFPED 

no 104 l8lfKK2 ■ ' ' 


DC 104 Ja3p4 
FDtJpI)**AtI»J5 


104 

106 


0 106 T»1 #FK2 
0 106 Js3i4 


DC 130 iJaliEK 


D a 0 

DC 120 I»lfK4 
IF (JCtUSCI) .NE 
DC 110 oa»lil0 
IF (NSTCJO) .NE, 
NCLUsiaai » J 


aj 

ij 


GO TO 120 
GO TO 110 


110 continue 
D a D + 1 
XltL) AtliS) 
X2C13 » AC1#4) 


12^ rCfJTINlJF 

rAT,L MIN(:n,L,FDta,6)) 
rAT.I. MAX(X2,L,Fn(J,5)) 

13P rnfjTiNUF 

WPITE (ft,3) kk 

no 132 as 1 , 4 
rnwsi ,PE+20 

TFfa,F.U,2,OP,a,EQ,45 OQWaO.fi 
no 1 32 1=1 ,KK 
132 Fr(I,a)*CON 

appangf and print nuTpiiT 
no 150 1=1, KK 
WRITE C6,4) I 
Waal ,/(K(6)+ENC3 • 

WP|TEC6,93 WJ 
WPITE(6,8) 

I a 0 

no 140 aal ,10 

IF (NfUlStai ,NE, I) GO TO 140 

L a L + I : 

MOLHCf ,l)=NAfJ) 

MASfli)sJ 

lF(MA(J3.F0.n 00 TO 136 
iF(n(j,n,LT,Fnci,n3 Fn(i,i)»B(a,n 
IF(B(a,2),OT,FI)CX,233 FD ( 1 , 2 ) bB CU , 2 ) 

GO TO 138 

136 ir(n(Jin,LT,FD(I,3)) FD ( T , 3 J=B (0 , U 
IF(R(J,23,GT.FO(I,4)3 FP ( I , 4 3 «B (U r 2 1 
138 cmniNUE 

140 CONTINUE , ■ . . , 

DO 142 JJ«1,D 
LI,«MAS(aJ3 

WRlTEt6,73 MCLU(I,JJ3,B(LL,53#B(Lt#7),TIMCLL3 
142 CONTINUE 
I1»I 

FDfII»ll«rDClI,13#WI13 
FDf Il,33»FDCII»33^tWf 13 
FDf JIr23aFDCn,23#WC23 
FD(TI,43»FDflI,43«Wf23 
FDin,53«FDCn»53«W(33 
Fn(II,63»FD(n»63«w(43 

WRITE (6, 6 3 rD(I,2 3,rD(I,n#FDtI|4 3#FD(I,3 3»FD(If63fFC(I,53 
150 CONTINUE , , 

200 CONTINUE 
STOP 
END 

"''subroutine BU11D(A»M,N,K,ITEB,SUM,NCLUS|DCLUS,X#NA,NN,MF,KK) 

dimension SUM(8,7,S5#A(10,75,XC73.NOLU5(103#DCLUSC103 
DIMENSION NAC10) 

DO 10 I«l»8 
DO 10 JaNNiN 
DO 10 KU*1,K 
10 SUMCIfafKU) « 0, 

DC IS KKK*1,KK 



w 


M 


74* 


SUMC2,a,KKK) 
St)Mr8,J,KKK7 s 1 
ei, 


= M^■ +KKK - 1 
nr IF* JSNF',N 

15 Sl'r(^»^»KKK) S A(KdK,vn 

nr f>0 Nr=t,TTFK 

I'iVftx s r, 

pr 10 KKK=1,KK 

nr ■^0 jxiSN,!'! 

ri0H(B,J,KKK) s 
IF (NC .EQ, 1) 

S!’HC2»»JfKKK) s , 

Sl'Md » JfKKK) s 

nr B0 ismh,h 
nr 40 JsNM,N 
X(J) = A(l,vn 

'^r A t L ^ K VE A N IR ( M » K K , SOM , X » NClIiS C 15 I r>CLl] S 1 1 ^ # Nin 

CONTINUE 
CCHTINUF 
rONTlNUf 
nr f*y kibukk 

Vf (fiUM( 2 ,;.Kn .NE. 0.5 GO to 76 

nr 75 i»tf8 
nr 75 J»NN,N 

SUMCI»<J,KI) = 0, 

75 CONTINUE 

76 CONTINUE 
140 CONTINUE 

return 
END 

subroutine KMEAN8(N,K,SUM»X,JN1N,DM1N|NN) 

DIHENSION 5UM(8#7,55 »X(7) 

JNIN ■ t 
PNIN » 

DO 20 a»i^K 
XP a l,0E«10 
DO • 0 » 

DO 10 l*NN»N 
XSUN * ABS tX (1 5 •• 

DD « DD ♦ XSUN«*2 
XP®XR4‘li 
CONTINUE 

no « SQPT(DD/XP5 
IF CDD »GT, DMIN) 

dnin » 00 

UNIN B U 
CONTINUE 

XN ■ N 

SUM£ if I»'3MIN5 * 'SUHOfliJ I 5 
40 CONTINUE _ ^ ■ 

return ' ,. . ' 

END -■ ' -ss* 


10 


SUMCl|IfU55 


GO TO 20 


20 


30 


n n 



c 



fUlHPOIlTTtHF SINGLEfX, COUNT, AVE,SD,XM1N,XM,AX»S5G) 

IV (COijKT ,PK, «,) GO TO 10 
AVB' = 0 , 
sr = 0 , 

Xf'lN = 1,0E + 20 
X^AX = - l, 0 E + l >0 
KPQ = 

If cril^T X COUNT + 1, 

AVF X aVE + CX-AVF3/CnilNT 
XAVF X ABSCX • AVF) 

ir (COUNT ,NF, 1,3 SSG a SSG COUNT « XAVR^nf 2/ (COUNT 
S^Q = SCJRTCSSO/COUNT) 

IF fXUIN ,GT, X3 XNIN .s X 
IF rXNAX ,LT, X3 XMAX a .x 
UFTUHN 


FNU 



5 


10 


5 


10 


SI'UPOUTINF N1N(X,NW,DMIN) 
nlMFNSlON X(10) 

TF CMM ,NF, 03 GO TO 5 
DNIN B 0 , 

BFTUPN 
CONTINUE 
PNIN » 1 ,0E + 20 
DO 30 
XO « X(«J3 

IF CXO ,GT. DMIN3 GO TO 10 

PVIM « XJ 

CONTINUE 

PFTUPN 

END 


SUHPOUTINE HAX(X,MM,DMAX) 
dimension X(103 
IF (MM ,NE, 03 GO TO 5 
DMAX a 0, 

BFTUPN 
CONTINUE 
DMAX * •! , 0 F + 20 
DO to J»l#MM 
XO » XCJ3 

IF (Xa ,LT. DMAX5 GO TO 10 

DMAX ■ %4 

CONTINUE 

RETURN 

END 
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